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Preface 

It  has  been  noted,  in  the  course  of  study  of  geodesy,  that 
there  is  a  spars ity  of  material  devoted  to  the  mechanics  of  determining 
the  size  and  shape  of  the  reference  ellipsoid  from  actual  observations. 
Bare  outlines  of  some  methods  are  given  in  the  few  current  textbooks 
available  in  English,  and  although  some  of  the  older  texts,  notably 
Clarke  (5)  and  Crandall  (6),  do  give  examples,  the  methods  used  are 
rather  ancient  and  outmoded  and  in  notation  systems  quite  unfamiliar 
to  the  present-day  student.  Hayford's  original  study  of  1909  is  the 
best  available  (10).  Since  one  of  the  basic  aims  of  geodesy  is  to 
determine  the  figure  of  the  earth,  it  seemed  like  a  good  idea  to 
pursue  the  matter,  if  even  yet  somewhat  superficially,  to  the  point 
where  there  would  be  few  questions  left  concerning  the  manner  of  going 
about  the  task. 

The  writer  wishes  to  acknowledge  the  help  and  advice  provided 
him  in  preparation  of  this  thesis,  particularly  by  his  adviser, 
Dr.  Wo  A.  Heiskanen,  Director  of  the  Institute  of  Geodesy,  Photo- 
grammetry,  and  Cartography,  the  Ohio  State  University.  Acknowledgement 
is  also  due  Dr.  U.A.  Uotila  and  Mr.  Ivan  Mueller,  both  of  the  Ohio 
State  University,  for  certain  helpful  suggestions;  the  U.S.  Army  Map 
Service  and  Mr.  W.M.  Kaula,  for  data  on  points  in  the  Finland  - 
South  Africa  arc|  the  U.S.  Navy,  without  the  instigation  of  which, 
the  course  of  study  referred  to  in  the  opening  sentence  would  not  have 
been  possible. 
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Finally,  the  writer  wishes  to  express  his  appreciation  to 
Frauke  and  Peter  Wilson,  the  former  for  her  fine  job  of  typing  the 
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Chapter  1.   Introduction 

1.1  Preleminary  Remarks. 

That  the  most  convenient  practical  reference  surface  for  purposes 
of  geodetic  computation  is  the  ellipsoid  of  revolution  has  long  been 
recognized.  This  is  emphasized  by  the  fact  that  the  major  datums  in 
use  in  the  world  today  are  all  based  on  ellipsoids  of  revolution.  The 
fact  that  this  shape  is  a  relatively  simple  regulai  mathematical  sur- 
face, and  that  it  closely  approximates  the  mean  geoidal  form,  make  its 
choice  quite  understandable.  The  Russians,  who  have  been  proponents  of 
a  triaxial  ellipsoid,  do  not  even  consider  its  use  for  practical 
geodetic  computations.   Izotov  (l8)  states s- 

"The  triaxial  ellipsoid  could  not  be  accepted  for  triangulation 
processing,  as  this  would  have  made  the  working  out  of  geodetic 
measurements  too  complicated". 

To  dismiss  the  subject  of  triaxiality  for  purposes  of  this  paper,  the 
following  is  cited  from  "The  Earth  and  its  Gravity  Field",  by  Heiskanen 
and  Vening  Meineszs- 

"It  appears  that  we  had  better  forget  the  triaxial  ellipsoid, 
particularly  because  since  the  earth  is  in  close  isostatic 
equilibrium,  it  can  hardly  have  the  shape  of  a  triaxial 
ellipsoid."  (13)  p.  79-80. 

The  use  of  geodetic  data  to  determine  relatively  precisely  what 
size  and  how  much  flattened  the  reference  ellipsoid  should  be,  is  the 
subject  of  this  discussion.   In  this  connection,  one  of  the  principal 

purposes  of  geodesy  has  been  variously  stated; - 
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"Its  theoretical  function  is  to  determine  the  size  and  shape 
of  the  earth  ..."  (13 )l  "Geodesy  is  the  science  which  treats 
of  investigations  of  the  form  and  dimensions  of  the  earth's 
surface."   (17)  J  "The  accurate  determination  of  the  figure  of 
the  earth."'  (6)|  "In  combination  with  observations  for  latitude, 
longitude,  and  gravity  to  assist  in  determining  the  size  and 
shape  of  the  earth,  ..."  (2);     "It  may  be  said  that  the  ultimate 
aim  of  scientific  geodesy  is  to  determine  the  size  and  shape  of 
the  geoid,  ..."  (2). 

Of  course,  data  collected  in  geodetic  surveys  along  with  the  calcula- 
tions connected  therewith,  and  computations  and  results  of  computations 
of  the  figure  of  the  earth  are  mutually  interdependent. 

Although  the  problem  of  determining  the  best  size  and  shape 
for  the  reference  ellipsoid  has  decreased  in  importance  relative  to 
other  modern  forms  of  geodetic  endeavor,  there  are  still  significant 
efforts  being  made  in  this  field,  notably  by  the  U.S.  Army  Map  Service. 
This  is  not  to  say  that  the  exact  dimensions  of  the  earth  are  not  im- 
portant, but  the  position  has  now  been  reached  where  further  precision 
given  to  a  reference  ellipsoid  must  be  subordinated  to  the  problem  of 
determining  the  deviations  of  the  actual  form  of  the  geoid  with  respect 
to  the  reference  surface  selected. 

This  change  of  emphasis  is  described  somewhat  succinctly  by 
Tengstrom; - 

"The  principal  problem  of  geodesy  is  to  determine  the  geocentric 
coordinates  of  all  points  on  the  earth's  surface,  a  task  which 
was  earlier  referred  to  as  a  determination  of  the  'size  and 
shape*  of  the  earth."  (30); 

or  by  L.  G.  Simmons  s- 

"The  ultimate  goal  of  the  geodesist  is,  I  suppose,  to  determine 
the  parameters  of  an  ellipsoid  of  revolution  which  best  fits  the 
figure  of  the  earth  as  a  whole.  But  this  is  not  all.  He  also 


is  concerned  with  the  details  of  lack  of  fit  of  this  ellipsoid 
to  the  actual  Earth's  shape,  the  geoid."  (28), 

Most  important  would  be  the  selection  of  one  reference  surface 
to  be  used  in  problems  of  world-wide  extent,  such  as  for  requirements 
of  missile  warfare,  the  success  of  which  might  very  well  be  determined 
by  accurate  knowledge  of  locations  of  launch  sites  and  targets.  Natu- 
rally, the  best  reference  surface  obtainable  is  to  be  preferred.  The 
establishment  of  a  "World  Geodetic  System"  and  conversion  of  existing 
systems  to  it,  has  for  several  years  been  one  of  several  significant 
projects  promoted  and  encouraged  by  Dr.  Heiskanen  (l6)  (13)  p.  299-310. 
Both  the  U.S.  Air  Force  and  the  U.S.  Army  Map  Service  are  developing 
world  geodetic  systems  for  military  purposes.  These  are  eventually  to 
be  combined  into  a  single  unified  Department  of  Defense  World  Datum, 
the  contributions  to  the  establishment  of  which  will  be  made  by  all 
U.S.  military  services.  The  reference  ellipsoid  to  be  used  for  this 
will  be  based  on  all  available  information  from  all  parts  of  the  world. 

(3)  p.  77. 

Having  exposed  the  near-obsolescent  character  of  the  problem 
of  determining  the  dimensions  of  the  reference  ellipsoid,  we  may  now 
proceed  to  describe  the  several  variations  of  its  solution. 

1.2  General  Discussion  of  Methods. 

The  standard  approach  to  the  problem  under  consideration  has 
been,  from  earliest  times  right  to  the  present  day,  that  contained  in 
the  arc  method.  Variations  of  this  approach  have  become  refined  and 
quite  sophisticated,  taking  into  account  the  earth* s  isostatic 


equilibrium,  and  in  some  cases  being  combined  with  the  gravimetric 
method.  The  earliest  recorded  application  of  the  arc  method  is  told 
in  the  now  trite  story  of  Eratosthenes*  activities  in  this  field. 

The  arc  method  assumes  that  short  meridional  arcs  on  the 
earth's  surface  are  arcs  of  circles;  when  two  meridian  arcs  are 
measured,  one  at  relatively  high  latitude,  and  one  near  the  equator, 
both  the  equatorial  radius  and  the  flattening  can  be  computed.  The 
more  such  arcs  used,  combined  in  a  least  squares  solution,  the  better 
accuracy  is  obtained. 

A  more  convenient  adaptation  of  the  arc  method  has  been  used 
since  early  in  the  present  century.  This  uses  a  combination  of 
geodetic  and  astronomic  data  in  the  form  of  astro -geodetic  deflections 
of  the  vertical.  Hayford  used  this  method  in  developing  what  is  now 
known  as  the  International  Ellipsoid  (10).  He  called  the  procedure  he 
used  the  "area  method",  because  the  data  was  taken  from  many  points  in 
a  large  network  of  triangulation.  He  made  several  solutions,  the 
deflections  of  the  vertical  in  each  being  reduced  to  the  geoid  (co- 
geoid)  according  to  the  theory  of  isostasy  using  different  depths  of 
compensation.  He  chose  the  solution  in  which  the  sum  of  the  squares 
of  the  residual  deflections  after  adjustment  was  the  least.  Thus  a 
"by-product"  of  his  work  was  a  solution  of  the  mean  depth  of  compen- 
sation for  the  earth's  crust. 

The  application  of  the  gravimetric  method  to  the  problem  under 
consideration  is  really  a  combination  of  the  arc  method,  modified  by 
the  introduction  of  gravimetrically  determined  deflections  of  the 


vertical,  to  obtain  the  earth's  equatorial  radius,  and  the  deter- 
mination of  the  earth's  flattening  from  gravity  anomalies  by  use  of 
Clairaut's  formula. 

The  latest,  and  of  course  the  most  modern  method  used  to 
determine  the  earth's  shape,  is  the  artificial  earth  satellite  method. 
By  measuring  the  perturbations  in  the  satellite's  orbit  around  the 
earth,  caused  by  deviations  from  the  spherical  in  the  earth's  shape, 
the  flattening  is  obtained.  ,  This  method  is  so  new  in  actual 
application,  that  an  evaluation  of  results  thus  far  obtained  can  only 
be  the  subject  of  conjecture. 

One  other  method,  similar  to  the  deflection  of  the  vertical 
application  of  the  arc  method,  involves  the  use  of  the  geoid  un- 
dulations, obtained  both  graviraetrically  and  astro-geodetically.  The 
solution  is  obtained  in  reducing  the  sum  of  the  squares  of  the 
differences  between  the  two  at  many  points  to  a  minimum. 


Chapter  2.  Dimensions  from  Two  Arcs 


2.1  Meridional  Arcs. 


In  this  analysis,  as  in  all  subsequent  ones  in  this  paper,  the 
reference  surface  is  assumed  to  be  an  ellipsoid  of  revolution.  The 
mathematical  properties  of  an  ellipse  form  the  basis  for  expressions 
referring  to  the  meridian  of  sucn  a  surface. 

Tne  measurement  of  two  arcs  in  meridian,  one  near  the  pole  and 
the  other  near  the  equator,  provides  the  data  for  this,  the  simplest 
means  whereby  the  earth  dimensions  may  be  obtained.  The  arcs,  which 
may  be  on  different  meridians,  are  measured,  and  the  difference  in 
latitudes  between  extremities  for  each  is  computed,  the  latitude  of 
one  end  point  for  each  arc  having  been  determined  astronomically. 


Meridian 


Fig.  1 
In  fig.  1,  cp.  and  cp  are  the  mid-latitudes,  and  s..  and  s  the 
measured  lengths  of  the  two  arcs.  Assuming  s.  and  s  to  be  arcs  of 
circles,  we  can  writes- 


sl 


(2-1)       R,  =  t—  esc  1" 
N    '         1    Acp, 

S2 

R0  =  -7-=-   CSC  1" 

2    Acp2 

Acp..  and  Acpp  being  expressed  in  seconds. 

The  expressions  for  the  radii  of  curvature  of  the  elliptical 
meridian  at  the  respective  latitudes  are: 

(2-2)        R     <*/•+*> 


1    ,,   2.2  sVz 
(1-e  sin  cp, ) 

R  =    a(l-e2) 


2  (l-e2sln\)S/* 

Equating  the  two  expressions  (2-1)  and  (2-2),  we  obtain: 

(2-%)  —   esc  1"     -     a(l'e  ) 

{d  5)  Acp.  CSC  J   ",,2.2  ,Vz 

\L  (1-e  sin  cp  ) 


S2     , h       a(l-e2) 
^2  (1-e  sm  cp2) 

Dividing  R.    by  Rp  in  the  form  (2-3)   and  solving  for  e    : 

,oM  ^2  l-e2sin2<p  2 

<2^  -^pT       "        ,2.2  =     q 

2     1  1-e  sm  cp 

2=    _i^2 


e     = 


2  2.2 

sm  9p-q  sin  cp 


Then,   the  value  of  a  can  be  found  from  (2-3) : 

sn  ( l-e2sin2co  f/z  bJ  l-e2sin2cp0)3/i 

(2-5)  a     =     -= - =     — - 

Acp1(l-e2)sin  1"  Acp2(l-e2)sin  1" 
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Now  all  that  remains  is  to  determine  the  flattening,  which  is 
simply: 

(2.6)    t,2±  =1.  _|_    .   i-v'— i- 

2.2  Oblique  Arcs. 

In  the  early  efforts  of  measuring  the  earth's  figure  relatively 
short  arcs  were  measured,  almost  always  in  meridian.   In  later  deter- 
minations, existing  triangulation  data  was  put  into  use,  and  consequent- 
ly successive  points  used  were  not  always  on  the  same  meridian.  Such 
geodetically  computed  arcs,  chosen  so  as  to  he  inclined  to  the  meridian 
at  as  small  an  angle  as  possible,  could  be  reduced  to  the  meridian  by 
appropriate  formulae  or  from  tables  based  on  the  used  ellipsoid.  Thus 
armed  with  the  same  data  as  in  section  2.1  a  similar  problem  can  be 
solved. 

When  the  length  of  the  oblique  arc  is  measured,  along  with  its 

azimuth  and  mid-latitude,  its  meridian  distance  is  given  by: 

12    2 
(2-7)     AcpE  sin  1"  =  -  s  cos  A  -  pTj—  S  sin  A  tan  cp 

1    S    2  2 

+  — ^  SJ  sin  A  cos  A  (l+3tan  cp) 

where  N  is  the  radius  of  curvature  of  the  prime  vertical  and  A  is  the 
azimuth,  reckoned  positive  clockwise  from  south. 

Using  (2-7)  a  chain  of  triang.es  can  be  projected  to  the  meridian 
for  the  whole  length  of  the  chain.  Given  geodetic  latitudes,  however, 
use  can  be  made  of  specially  prepared  table,  such  as  (3^)>  (35),  or 
(36),  which  have  meridian  distance  tabulated  against  latitude  for  the 
ellipsoid  involved. 


The  above  method  of  determining  the  earth's  dimensions  (Sections 
2.1  and  2,2),  though  it  represents  the  basis  for  other  methods,  is 
mainly  of  historical  and  mathematical  interest,  since  no  longer  are 
serious  scientific  computations  carried  out  in  such  simple  terms. 
References  for  the  presentation  given  in  these  two  sections  are  (6), 
(17),  and  (23). 


Chapter  3«   Dimensions  from  Several  Arcs. 

3.1  General. 

When  geodesists  of  the  l8th  Century  first  began  to  compare  the 
results  of  their  arc  measuring  efforts,  they  found  such  large  discrep- 
ancies that  they  began  to  speculate  about  the  accuracy  of  their  meas- 
urements and  about  the  correctness  of  the  assumption  that  the  earth's 
true  shape  is  an  ellipsoid  of  revolution.   Of  course  the  accuracy  of 
surveys  of  those  early  days  was  not  what  it  is  today.   Nor  could  one 
hope  to  obtain  the  same  results  from  arcs,  taken  two  by  two  in  dif- 
ferent parts  of  the  world  having  varying  kinds  of  topography,  meas- 
ured by  the  crude  methods  of  the  day.   Not  only  was  the  measurement  of 
geodetic  triangles  so  rude  that  spherical  excess  remained  undetected, 
but  the  process  of  adjustment  by  the  method  of  least  squares  was  un- 
known (23 )  • 

Obviously,  not  wanting  to  discard  such  a  convenient  reference 
surface  as  an  e'llipsoid,  which  after  all  departs  from  the  geoid  sur- 
face by  probably  less  than  100  m.  at  the  maximum,  and  from  the 
equilibrium  figure  by  even  less  (of  the  order  of  2  m.  at  cp  =  0,^5  > 
and  90,  (13)  P-  59)  ?  "the  best  thing  to  do  would  be  to  combine  the 
results  of  several  arc -measurements  in  such  a  way  as  to  determine  the 
mean  values  for  the  earth  as  a  whole. 

3-2  Laplace's  Method. 

Toward  the  end  of  the  l8th  Century,  Laplace  published  the  re- 
sults of  his  attempt  at  combining  the  data  from  nine  different  sources, 
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extending  from  Lapland  to  Cape  of  Good  Hope  in  latitude.  His  method 

involved  the  data  corresponding  to  an  arc  of  one  degree  of  latitude  for 

each  of  nine  different  meridians. 

Again  considering  these  one  degree  arcs  as  arcs  of  circles  and 

using  R,  the  radius  of  curvature  of  the  midpoint  of  the  arc,  and  cp, 

the  mid-latitude,  we  havet- 

/o  ,  v      -   2  Ti  R        7T  a  (l-e  ) 
(3_1)     *  "  ^°"  =  !80  (l-e2  sin2 >* 
Developing  this  in  series,  we  obtains - 

/,  0n      ,   IT  a(l-e2)   ,..   3  2   .  2     15  k     .   k  s 

(3-2)     d  =  - jg^ — '-     (1  +  |  e  sm  cp  +  ^  e  sin  cp  +  . . .  ) 

which  can  be  written  in  the  following  form:- 

2         k 
(3-3)     d  =  M  +  V.   sin  cp  +  P  sin  cp  +  . . . 

j   uj  u        7ra-  (l-e<~) 

in  Which  M  =   yTTTr — - 

3     2 
N  =  |  e     M 

15     k 
P  =  rp  e     M,    . . .    etc. 

Powers  of  e  higher  than  the  square  are  neglected. 

In  order  to  combine  the  data  for  each  degree -measurement  in 
the  manner  yielding  the  least  discrepancies,  the  observation  equations 
for  each  arc  were  of  the  forms - 

(3-k)  d  -  M  -  N  sin2  cp  =  v 

2 
dp  -  M  -  N  sin  cp  =  v?     etc. 

in  which  v  ,  v  ,  etc.,  are  the  residual  errors.  Since  the  method  of 
least  squares  was  unknown  at  the  time,  Laplace  combined  the  error 
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equations  in  such  a  way  that  the  absolute  values  of  the  residual  errors 

would  be  a  minimum,  and  the  algebraic  sum  equal  to  zero.  Then  solving 

2 
simultaneous  equations  for  M  and  N,  the  values  of  e  and  a  are  found. 

Of  course,  this  problem  is  easily  adapted  to  least  squares  solution. 

The  results  of  Laplace's  computation  ares- 

M    =  11  06l.2  meters 
N     =     1  19^.9  meters 
e2  =  0.007202 
a     =     6  383  600  meters 
f     =  1/278 

The  error,  or  difference  between  a  measured  degree  and  one 
computed  from  the  result  of  the  crude  adjustment  of  the  nine  arcs,  is 
of  the  order  of  268  meters.  From  this,  Laplace  concluded  that  the 
earth  deviates  materially  from  an  elliptical  figure  (23). 

3<>3  By  Correction  of  Latitude. 

A  refinement  of  Laplace's  method,  in  use  by  the  end  of  the  19th 
Century,  is  here  described  (6),  (17).  This  also  involves  the  formation 
of  observation  equations  between  observed  latitudes  and  either  directly 
measured  meridian  arcs  or  oblique  arcs  projected  to  the  meridian. 

For  short  arcs  the  relation  s  =  Acp  R  sin  1"  can  be  used,  but 
for  longer  arcs  a  correction  is  required.  Crandall  (6)  states  that  the 
simple  relation  given  is  good  for  a  Acp  of  "several  degrees"  because  of 
the  large  probable  error  of  latitude  determination,  some  OVO^  (more 
likely  this  will  be  on  the  order  of  O'.'l  at  best),,  in  addition  to  the 
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error  caused  by  local  deflections  of  the  vertical.  The  correction,  for 
use  with  longer  arcs  is:- 

(3-5)  d  s  =  g  a  e2   (Acpst  sin  l")3  cos  2  cp 

Given  astronomic  latitudes  for  several  stations  in  different 
arcs,  as  well  as  geodetically  measured  segment  lengths  s_,  we  choose 
separate  initial  points  for  each  arc.  From  the  first  two  latitudes  in 
any  one  arc  we  have:- 

(3-6)     cp2  -  cp1  =  p— ^  -  R  s?nV 

where 

/-.  2.2      N3/2 

(3-7)  |  =  (l"e     sin2^ 

R  a  (1  -  e2) 

In  this  method  we  make  use  of  the  used  reference  ellipsoid,  for 

2 
which  we  know  the  values  a  and  e  ,  and  compute  corrections  5a  and 

2 

8e  ,  to  that  ellipsoid  resulting  from  the  data  given.  Thus  a  =  a  +  6a, 

2     2     2 
and  e  =  e   +  oe  . 
o 

Expanding  equation  (3-7)  into  a  series  and  substituting  it  in 
equation  (3-5) J- 


<3-8)    %  "  «*-  iET 


1         5a  ,.        3     ,   2    v   6e2 

s 2 +  &  -  § sin  *>  — 

o       a  o 

o 


R     sin  1" 


2 
in  which  R  is  the  value  of  R  for  a  and  e  . 
o  oo 

Considering  the  measured  meridional  arcs  accurate  in  comparison 
with  the  observed  latitudes,  we  give  corrections  v.,  v  ,  to  latitudes 

cp.,  cp  ,  etc.  Since  in  the  derivation  of  (3-8)*  all  terms  containing  e~ 

2        ? 
were  neglected  in  the  coefficients  of  5  a  and  8  e  ,  i.e.,  e  was  con- 
sidered equal  to  zero,  we  can  say  R  =  a  and  s/(a  sin  1")  =  cp  -  cp  . 


Ik 


Thus  equation  (3-8)  becomes;- 

(vv2)  -  <vV  -  htht  +  rnrr  [- 15  ^  i-i»W] 


S 

[- 

r  i 

8a 
a 

0 

and 

then 

in  i 

a  i 
o 

3in  1" 
d  s 

R 
o 

sin  1* 

o 


9o  "  9n  a  2 


(3-9)  v2  =  Tj  -  a  2sin  p  5  a  +  (1  -  |  sin2  cp)  (cp2  -  q^)  5 


e 


s  -  ds    /       \ 

+  R  sin  1"  "  <«fe  "  V- 
o 

This  becomes  the  basic  observation  equation  to  be  written  for 

2 
each  segment  of  each  arc.  The  unknowns  are  v^ ,  o  a,  and  5  e  . 

2 
In  the  least  squares  solution,  in  which  Z  v  =  minimum,  the 

magnitudes  of  the  various  coefficient  elements  make  it  convenient  to 

use  the  following  modifications. 

Placing  x  =  5  a/ 1000       y  =  1000  8  e 


nu  = 


=  -  1000  (cp2  -  91)/aQ 

n2   1600-  (1  -  2  SlD  CP) 

i      s  -  ds     /       s 

\   =  R  sin  1"   "  (cp2  *  V 
o 

We  can  then  write  a  simple  error  equation  for  each  segment, 
substituting  the  above  symbols  in  (3-9)?- 

v2  =  vx   +  m2  x  +  n2  y  +  1^ 

Since  the  initial  latitude  must  be  corrected,  an  additional 
observation  equation  is  written  for  it,  simply  V-  =  v. .  In  this 
method,  there  will  be  as  many  unknowns  as  there  are  arcs  plus  two,  i.e., 
the  v  for  the  initial  latitude  of  each  arc,  plus  x  and  y.  In  the  case 
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that  a  value  for  the  initial  v  is  assumed  to  be  some  arbitrary  value, 
e.g.  zero,  there  will  be  only  two  unknowns 0 

For  each  arc,  observation  equations  will  be  as  follows  - 

(3-10)  v1  =  v± 

v1  +  m2x  +  n2y  +  lg  =  vg 

v..   +  m„x  +  n_y  +  1_  =  v_ 
1         j  j  5         j 

*  *  #  *  # 

vn   +  m  x  +  n*y  +1     =  v 
In  ir         n         n 

The  result  of  a  computation  of  this  type  presented  by  Crandall 

(6),  and  using  data  (circa  1890)  given  in  the  "Handbuch  der  Vermessungs- 

kunde"  (37)  for  six  arcs  (French,  English,  Hanover,  Prussian,  Russian, 

and  Swedish),  varying  in  length  from  0.5  to  12.4,  is  as  follows ;- 

a     =  6  377  800  meters 
e2  -  O.OO6909I 
f     =  I/289 

The  standard  error  of  the  mean  amounts  to  2'.'1  for  a  latitude 
determination,  which  is  quite  a  bit  greater  than  the  accuracy  obtain- 
able for  first  order  work,  (about  O'.'l).  Thus,  Crandall  concludes  (6), 
an  ellipsoid  of  revolution  will  not  fit  the  data  without  large  local 
deviations  of  the  vertical. 

3.4  By  Correction  of  Radius. 

A  method  given  by  Bomford  (2)  p.  3^4,  which  is  a  variation  of 
that  described  in  Section  3»3>  works  with  radii  of  curvature  at  points 
in  an  arc  rather  than  with  latitudes. 
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For  given  values  of  R  in  widely  separated  latitudes,  the  values 

2  /    x 

of  a  and  e  are  obtained  from  the  expansion  of  equation  (2-2 )t~ 

(3-11)   R  =  — —  o  — p~ — Tfo 
(1  -  e2  sin2  cpr7 

2       3  2   3  4    4 
=  a  (l  -  e  )  (l  +  ^  e  +  tt  e  sin  cq  +  . . . ) 

2         2  2 

Using  e  =  2  f  -  f  and  the  approximation  6  e  «  2  6  f,  and 

2 
dropping  the  terms  of  higher  order  than  >  ,  the  following  form  is 

obtained i - 

(3-12)        5  R  =  6  a  -  2a  5  f  (l  -  J  sin2cp) 

O  c. 

in  which  8  R  =  R  -  R  results  in  changes  5  a  and  6  f  in  a  and  f   1 
m    o  ^^  o     o 

the  ellipsoid  provisionally  adopted.  From  equation  (2-3)  we  may  use 

R  =  s  esc  r'/^p"  for  arcs  shorter  than  about  100  miles  (2).  The 
m 

value  of  R  would  be  obtained  for  the  used  ellipsoid. 

In  a  more  practical  application,  more  than  two  values  of  R 
would  be  available,  and  observation  equations  for  least  squares 
solution  would  be  as  follows- 

S  CSC  l" 

(3-13)  w±  ,  »  .  2aoy  (1  -  |  sin2  ^)  -  ^H   ♦  (Ro)1 

s  csc  1" 
v£  .  x  -  2aoy  (1  -  |  sin2  cp2)  -  -^r-  *   (RQ)2 

#     #    *    #    # 

op  s    CSC    1" 

v     =  x  -  2a  y  (1  -  4  sin     9   )   "     "A      » —  +   (R  ) 
n  oJ    v         2  Yny  Acp  v  o'n 

Yn 

where  x  =  5  a  and  y  =  6  f  are  the  unknowns,   s     is  the  measured  length 

p 
of  the  arc  P«Pp5   and  Acp  "  =  cp     -  cp_    in  seconds,,   etc.     Again  Z  v     is 

made  a  minimum.  Weights  may  be  assigned  as  deemed  necessary. 


Chapter  k.      Dimensions  from  Arcs,  Using  Deflections 

of  the  Vertical 

4.1  Basic  Considerations  -  Astro-Geodetic  Deflections. 

(Note;  A  more  complete  analysis  of  "basic  considerations,"  in- 
volving absolute  deflections  of  the  vertical,  is  given  in  section  ^-.3») 

As  mentioned  earlier,  a  very  convenient  adaptation  of  the  arc 
method  makes  use  of  astro-geodetic  deflections  of  the  vertical,  which 
are  generally  available  in  profusion  wherever  triangulation  control  nets 
nave  been  made.  This  has  been  done  in  many  parts  of  the  world  in 
connection  with  the  various  national  surveys  being  carried  out. 

This  method  can  be  applied  to  individual  arcs  of  meridians  in 
one  extreme,  and  to  the  whole  network  of  triangulation  points  in  the 
area  covered  by  a  national  control  survey  in  the  other.  Because  of  the 
latter  application  of  the  method,  Hayford  chose  to  call  it  the  "area 
method"  when  the  results  of  the  U.S.C.  &  G.S.  computations  of  the  figure 
of  the  earth  based  on  data  collected  in  the  large  area  covered  by  the 
United  States  were  published  (10). 

Here  again,  the  reference  ellipsoid  used  in  computation  of  the 
geodetic  data  is  used  as  a  point  of  departure,  or  more  accurately,  an 
"area  of  departure",  from  which  we  must  apply  measurements  and  compu- 
tations to  get  as  close  to  the  best  ellipsoid  as  we  can.  This  "best" 
ellipsoid  is  the  one  most  nearly  approximating  the  earth  spheroid  or 
mean  geoid  surface,  in  size,  shape,  and  orientation. 

Figure  2  shows  in  simple  terms,  the  relationship  which  can  exist 
between  a  reference  ellipsoid  and  the  geoid. 
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Fig.  2  -  Astro-Geodetic  Deflections  of  the  Vertical 

In  general,  there  are  six  different  items  that  may  he  corrected 
by  the  procedure  involved,  and  the  particular  ones  that  are  corrected 
depend  upon  what  is  desired  and  upon  what  and  how  much  data  is  supplied. 
The  existing  geodetic  data,  i.e.,  latitudes,  longitudes,  deflections, 
etc., depend  upon  the  initial  point  (datum)  of  the  survey,  and  upon  the 
assumptions  made  about  that  point.  The  six  items  are:  the  equatorial 
radius  a;  the  polar  flattening  f ;  and  the  following  quantities  assumed 
for  the  initial  point:  its  latitude  cp  ,  its  longitude  \   ,  the  geoid 
height  M  ,  and  the  azimuth  A  of  its  meridian.   In  many  cases,  all  that 
are  desired  would  be  the  corrections  £>a  and  Sf.  Usually  the  cor- 
rections to  cp  and  X  at  the  initial  point  will  be  in  the  form  of 
Yo      o 


corrections  to  the  meridian  and  prime  vertical  components  to  the 
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deflections  of  the  vertical. 

The  following  notation  will  he  used  in  the  discussion  to  follow: 
cp'  ,   \'  -  Astronomic  latitude  aad  longitude 

cpj  X  -  Geodetic  latitude  and  longitude 

(cp  -cp)     -  Meridian  component,  astro-geod.  defl. 

(A.  -\)coscp  -  Prime  vertical  comp.,  astro-geod.  defl. 


f 
v 

dVD 

dio 


-  Residual  deflection  in  meridian 

-  Residual  deflection  in  prime  vertical 

-  Orientation  corrections  to  deflections 

due  to  displacement  of  ellipsoid 

-  Change  of  ellipsoid  corrections  to  deflec- 
tions. These  contain  the  values  Sa 
and   Sf  which  result  when  an  ellipsoid 
is  changed.  These  corrections  actually 
apply  to  the  geodetic  parts  of  the  astro- 
geodetic  deflections. 


When  it  is  desired  to  find  a  better  orientation  of  the  ellipsoid 
as  well  as  to  find  its  size  and  shape,  the  observation  equations  which 
may  be  written  for  each  point  take  the  following  general  form: 


(4-1) 


£     =     (cp'  -  cp)  +  d 


D  +  dfC 


>1      =   ( \J    -  X)    cos  cp  +  d^D  +  d^/c 


If  it  is  desired  only  to  find  a  and  f ,  using  data  from  a  meridian 
arc  only,  the  series  of  observation  equations  will  be  as  follows: 
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(4-2)         i±  =  (cPl'  -  cp1)+  (d  ic)1 

i2  =  (<P2'    -  ro2)+  (d  lc)2 
#  *  *  *  * 

|      =  (cp  /    -  cp  )   +  (d  |    ) 

2 
Then  the  least  squares  solution  will  involve  making  Z|  = 

minimum,  this  giving  us  the  parameters  a  and  f  of  an  ellipsoid  which 

will  best  fit  the  condition  which,  in  effect,  forces  the  residual  de- 

flections  to  a  minimum,  as  close  to  zero  as  possible.  Weights    '  be 

assigned  to  the  observation  equations  as  deemed  necessary. 

In  general  form,  the  condition  can  be  stated: 

v   2     2 
(^-3)    L(l   +  7  )  =  minimum, 

which  is  good  for  any  problem.  Either  or  both  terms,  £  and  >7  ,  can 

be  used  according  to  the  problem  to  be  solved. 

In  this  case,  in  which  astro-geodetic  deflections  are  used,  the 

better  results  are  obtained  when  these  deflections  have  been  reduced 

isostatically. 

k.2     Formulas  for  Displacement  and  Change  of  Ellipsoid  Corrections. 

There  are  four  well-known  methods  of  computation  of  the  change 
of  ellipsoid  and  displacement  corrections  to  the  astro-geodetic 
deflections.  These  are  based  on  formulas  representing  different- 
derivation  approaches  to  the  problem.  They  are  the  formulas  of  Hayford 
(10),  the  U.  S.  Army  Map  Service  (k) ,   Bomford  (2),  and  Vening  Meinesz 
( 33) •  The  derivations  of  these  formulas  are  given  in  the  references 
cited  (except  (^)),  and  being  beyond  the  scope  of  this  paper,  will  not 
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"be  covered  here. 

4.21.  Hay-ford* s  Formulas.  The  form  of  observation  equation  given 

in  (10)  is  as  follows: 

(k-k)        |  =  kx  (q>)  +  1±   (\)  +  m±   (A)  +  n^jjfc)   +  o1  (10,000  e2)  +  ( q. '  -cp) 

^  =  k2  (cp)  +  12  (X)  +  m2  (A)  +  n2^)  +  o2  (10,000  e2) 

+  (\f    -  \)    cos  cp 

in  which  (cp),  (X),  (A)  are  unknown  corrections  to  the  initial  cp,X,  and 

a  2 

A  at  the  origin  of  the  survey,  and  (^7577)  an<^  (10,000  e  )  are  unknown 

2 

corrections  to  the  values  of  a  and  e     for  the  reference  ellipsoid.      It 

can  he  seen  that  these  are  similar  in  form  to  equations  (4-1). 
The  change  of  ellipsoid  corrections  are: 

(1J.-5)      d£c   =  n1   (jgy    +  0l  (10,000  e2) 
d\  =  n2   (I50}    +  °2  (l°'°00  &2) 


100  N 

a  sin  1"    "-  R 


in  which  n_    =  — : — _ _.  _    .,  » — .— —  e     cos  A^ 


•  ~l 


2       2 
e      0  2  2 

1  +  75--  cos  <P     cos  A-, 

6(l-e2)  °  T 


.    2 

sin  op 

0     =  — ^ P Q   cos  A^ 

20,000  sin  1"   (l-e  sin  q>  ) 

-  -§  i *-«  (cp-cp J  [1+2  sin2cp     +  3  cos(cp  +cp)  ] 

H     40,000   sin  1"   (1-e^r  000 

„  .  _1?0  9  sin  A 

2  a  sin  1  a 

.    2 
sin  cp 

o     = 2 _ _ e    sin  Ag 

20,000   sin  l"(l-e  sin  cp  ) 
/  o 

(l-e2sin2cp  )1/2  cos2cpo  2  2 

+  5-= —     S  6      cos     L  sin  A^ 

60,000  a  sin  1"   (l-edr 


22 


The  corrections  for  displacement  (orientation)  are:- 

(k-6)  d|D  =  kx  (cp)  +  lx  (x)  +  i^  (A) 

dT)D  =  k2  (cp)  +  12  (X)  +  ra2  (A) 

c          sin  (X-X  ) 
in  which     k   =  -  1  +  -  (l  +  Q)  g-y 


2  sin  £CW 


1   =  zero 


c         sin  A^  [1  +  cos  (X  -  X  )] 

i  =  I  (  i +  Q  )  — ^-n — 

1    R  ^  sin2  |  (A3  -  V 


k   =  -  sin  cp  sin  ,(X  -  X  ) 
lp  =  -  cos  cp 


m. 


cos  cp  cos  A^  sin  (X  -  X  ) 


2  sin 


*F 


2 
Q   =  f-  cos2  |  (^  -  Ap) 

The  purpose  of  the  coefficients  k,  1,  m,  n,  and  o  can  be  stated 
as  in  this  example  for  one  of  them:  k..  is  a  numerical  coefficient, 
such  that  if  the  latitude  of  the  initial  point  were  corrected  by  the 
amount  (cp),  the  change  produced  in  cp1  -  cp  would  be  k.  (cp). 

It  can  be  seen  from  inspection  of  the  coefficients  that  for 
points  on  the  same  meridian  as  the  initial  point,  the  values  are 
k.  =  -1,  1  =  0,  and  m.  =  0,  so  that  the  observation  equation  would  be:- 

(h-1)  g  =  (cp)  +  nx  (^)  +  0l  (10,000  e2)  +  (q>«  -  cp). 

4.22  Army  Map  Service  Formulas.  In  an  article  in  Transactions  of  the 
American  Geophysical  Union,  No.  37,  1956,  Chovitz  and  Fischer 
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reported  on  the  U.  S.  Army  Map  Service's  determination  of  the  figure  of 

the  earth  from  arcs  ('  ).  The  procedure  is  described  in  some  derail, 

the  essential  part  regarding  correction  formulas  being  reproduced  in 

somewhat  different  form  here.  The  method  is  designed  to  be  correct  for 

meridional  arc  segments  of  1  ,  but  not  more  than  3  in  length,  and  for 

arcs  of  parallel  of  comparable  length. 

For  meridian  arcs  the  formulas  are: 

<M)   aec  .  [-|5(,n  -  „  )  ♦  6e2  {«„-,)  .£»  (%  .  „  ) 

L  n=l 

II  cos  2cp     -  III  cos  4cp     +  IV  cos  6cp     }     esc  1" 

L        m  Ym  YmJJ 

in  which  cp  =  mid-latitude  of  segment  =  —  (cp  _  +  cp  ) 
Ym  °      2   rn-l   Mr 

(4-8)   d£_  =  i     -(cp  '    -   cp  )  =  initial  orientation  of  the  arc  segment. 

In  this  method  the  computation  proceeds  segment  by  segment,  with 
a  new  initial  point  for  each  segment  of  the  arc  when  the  length  limits 
given  above  are  reached.  Thus  the  procedure  is  similar  to  the  one 
described  in  section  3»3«  This  presupposes  that  all  the  points  in  the 
arc.  furnish  data  that  is  adjusted  to  the  same  datum.  In  contrast,  the 
Hayford  formulas  (and  those  to  be  given  subsequently)  are  designed  to 
.over  a  large,  but  limited  area,  using  only  one  initial  point  in  the 

computation.  The  unknowns  in  formulas  (4-7)  and  (4-8)  are  |  ,  Sa,  and 

p 

Se   .  The  Roman  numeral  coefficients  are: 


IV  -  h'k 


II 
III  = 


1 
4 

7  2  17  4 
TB  e          32 

3 

IT 

3  2  ,  45  4 
±6             32 

2k 


For  arcs  of  parallel  the  formulas  are: 

zi^nj-^          ("Sa       -.-vn                    r  2,.      ,    N  sin     cp  cos  cp 

(4-9)      &yn  =     — —  ,\  -X  )    cos  cp  +  de  IX.  -X  J   — — 

v  "          '  C   L  a   *  n  o'            Y       no'  - 


/-,   2  .  2     k   ,  k   x 

^1+e  sin  cp  +  e  sin  cp) 


esc  1" 


where     cp  =  latitude  of  the  parallel. 

(4-10)  d^=Y"(^  ~  ^  )  coscp  =  initial  orientation  of  the  arc  segment 

2 

The  unknowns  are  ^  ,  £a,  and  Se  .   Solutions  for  meridians  and 

parallels  can  be  done  separately  or  can  "be  combined  for  a  cross -solution 

M. 

4.23.  Bomford's  Formulas.  The  formulas  for  change  of  ellipsoid 
given  in  Bomford  (2)  p.  127-130,  actually  contain  corrections  for  both 
displacement  and  change  of  ellipsoid.  They  are  developed  here  in 
separated  form.  These  formulas  are  used  In  Bomford's  version  of  the 
area  method,  which  supposes  a  "homogeneously  computed  national  survey" 
covering  a  20  by  20  area  containing  50  astro -geodetic  stations. 
Bomford  states  that  it  will  generally  be  best  to  ignore  values  of  >? 
deduced  from  azimuth  observations  (these  were  omitted  in  the  description 
of  Hayford's  method  above),  and  to  use  only  those  based  on  longitude 
(2)  p.  365. 

The  change  of  ellipsoid  formulas  as  given  by  Bomford  (2)  p.  129, 
365,  are  as  follows: 

(4-11)  d|  =-   'l/a)sir.  v.   cos  co  (U  sin  u  +  \  cos  u  ) 

+  sin  u  sin  to  bri     -  (l/a)cos  u(V  sin  u  -  U  cos  u  ) 
•  o     '  o        o 

+  sin  2u  5f  I  esc  1" 
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(4-12)   d>?  =  -  (l/a)  sin  cd  (U  sin  u  +  V  cos  u  )  -  cos  cd  S>j    1  cbc  1" 

where       U  =  -a  <5  |  -  a  <5  f  sin  2  u 

*o  o 

r      2 
V  =  N  -Sa  +  adf  sin  u 
o  o 

u  =  reduced  latitude;  tan  u  =  (l  -  f)  tan  cp 

eo  sb  longitude  measured  east  from  origin 

N  =  height  of  origin  in  corrected  ellipsoid 

above  that  in  reference  ellipsoid 

<5 i     and  b~n    are  in  radians, 
o  'o 

Substituting  the  values  of  U  and  V  in  (4-11)   and  (4-12),   and 

separating  the  expressions  into  two  parts:      d|  =  d|p  +  d£ 

d>7    =  d>?^+d>f,we  ohtain: 


(4-13)  d|c  = 


(sin  u  cos  cd  cos  u  -  cos  u  sin  u  )  5 a/a 
v  o  o'         ' 


.   2 
cos  a>  sin  u  cos  u  +  cos  u  sin  u 
o     o  o 


sin  u 
(l  +  cos  u  )  -  sin  2u  of 


esc  1" 


(4-14)  &£^   =  (sin  u  cos  ao  sin  u  +  cos  u  cos  u  )  <5  6  " 


sin  u  sin  a>  8)?"-   (sin  u  cos  a>  cos  u 
«o  v  c 


-  cos  u  sin  u  )(N  /a)  esc  1" 

o   o' 

(4-15)  d  ^p  =  sin  oo  cos  u  (  5  a/a)  +  sin  oo  sin  u  cos  u   Sf  esc  1" 
(4-l6)  d>7   =  sin  oo  sin  u  <5£"  +  cos  cd  S^" 

-  sin  cd  cos  u  (N  /a)  esc  1" 

o  o'  ' 

4.24,  Vening  Meinesz'  Formulas.  The  formulas  here  referred 
to  are  those  published  in  1950  (33)>  and  which  were  developed  for  the 
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purposes  of  (a)  describing  the  relationship  between  the  variation  of 
the  vertical  deflection  components  at  the  origin  of  a  survey  and  the 
changes  resu]ting  in  the  deflections  at  any  arbitrary  point  in  the 
survey,  the  variation  having  been  caused  by  a  parallel  displacement 
of  the  reference  ellipsoid  with  respect  to  the  geoid;  and  (b) 
describing  the  changes  in  the  deflections  of  the  vertical  brought 
about  by  a  change  of  reference  ellipsoid. 

The  complete  formulas  referred  to  in  (a)  above,  are  as 
follow: - 

(k-11)   d£  =  -R°  [cos(cp-cpo)-2  sin  cp  sin  crc  sin"  |(\-\  )]  6g 

r 

+  =-  sin  cp  sin  (\-X  )  5^1 
t\  oo 

1  2  1 

+  =•  [sin  (cp-cp  )  -  2  sin  cp  cos  cp  sin  ^  fX-\   )]&N 
R      VY  Yo'  o      2     oJ       o 

R  r 

(U-18)  dTi  = sin  cp  sin  (  \-\   )&|  +  —  cos  (\-\  )   5^ 

v    '   D    r     ^o     v    oJ   bo   r      v   oy   o 

cos  cp  sin  (  \-\   )  5N 

r     Yo     v    o;   o 

where  R  ,  R  =  radius  of  curvature  in  meridian  at  initial  point  P 
o  r  o 

and  at  point  P,  respectively 

r   r  =  radius  of  curvature  in  prime  vertical  at  initial 
o, 

point  P  and  at  point  P,  respectively  (r  is  used 
here,  rather  than  the  usual  N,  which  is  reserved 
for  geoid  distance). 

The  above  form  is  that  given  by  Heiskanen  in  (13)  p.  302,  and 

is  better  suited  to  the  use  of  tables  than  the  form  in  the  original 

2   •? 

paper  (33 )>  which  uses  a(l-e  )/w  rather  than  R,  and  a/w  rather  than 
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2   2  2- 
r,  where  W  =  (l-e  sin  cp)  . 

Vening  Meinesz  points  out  that  his  formulas  are  similar  to 
those  developed  by  Helmert  (31),  hut  include  on  terms  as  essential 
differences. 

If  the  line  P  P  is  long,  the  complete  formulas  above  must  be 
used.   However,  for  nets  of  small  east-west  extent,  terms 
proportional  to  the  squares  of  (\-\  )  can  be  neglected;  and  for 
nets  of  small  extent  in  latitude,  the  ratios  R^R*  r  /R,  R  /r5  and 
r  /r  are  very  nearly  1  (within  0.01),  and  the  formulas  (4-17)  and 
(4-18)  become s- 


(4-19)  d^  =  cos  (cp-cp  )8f  +  sin  cp  sin  (\-\  )§t\ 

8N 

+  sin  (cp-cp  )  — 
o  a 


(4-20)  dT)^  =  -  sin  cp  sin  (\-\   )  8g  +  Br] 
v    '      'D        Yo     x  o'     bo    'o 

5N 

-  cos  cp  sin  (\-\  )  

Yo     N   o'  a 


According  to  (13 ),   the  formulas  (4-19)  and  (4-20)  apply  when 

P  P<  300  km.  in  east-west  direction  and  <  1500  km.  in  north -south 
o 

direction.  However,  Vening  Meinesz  (33)  states  that  if  it  is 
desired  not  to  neglect  quantities  greater  than  l/300  of  the 
deflection  of  the  vertical,  these  formulas  can  be  used  for  nets 
extending  up  to  4  40'  in  east-west  direction  and  to  more  than  20  in 
latitude,  from  the  origin  P  .  Furthermore,  the  complete  formulas 
will  seldom  require  that  the  ratios  between  R,  R  ,  r,  and  r  be  any- 
thing but  1  even  for  large  nets,  to  attain  the  required  accuracy. 


28 

The  complete  formulas  referred  to  in  (b)  above,  for  change  of 
ellipsoid,  are  here  transcribed:;  - 

(4-21)  d|  =  {[sin  (cp-cp  )  -  2  cos  cp  sin  cp  sin2  |  (x-X   )]  66 

1  1 

-  4  cos  cp  cos  -  (cp+cp   )   sin  ^   (cp-cp   )  5f 

o 
+   [(2-3  sin     cp  -  cos   (cp+cp   ))  sin  (cp-cp   ) 

2  2  2  1 

-  2  cos  cp     sin  cp  (2+sin  cp     -  .3  sin  co)sin     =■  (X-X   ) 

tq  v   \  ^o  '  '  2  o' 

o  i  -i 

-  sin  2cp     sin     ^  (cp-cp   )   -  2  cos  cp  cos  ^   (cp+cp   ) 

sin  |  (cp-cpQ)]   f5B  +  [  |  sin  4  cpQ  (sin2  |  (cp-cpQ) 

2  1  3 

+  sin  cpo  sin  cp  sin     ^  (*--*-0))   -  slr>  cpQ  sin  -  (cp+cpQ) 

3  11 

sin  |  (cp-cpQ)   -  sin  cpQ  sin  -  (cp+cpQ)  sin  -  (cp-cpQ) 

1  1 

-  4  cos  cp  cos  -  (cp+cpQ)   sin  ^  (cp-cpQ) 

2  1  2  1  T 

+  4  sin  2  cp  cos     -  (cp+cp   )  sin    '  -  (cp-cp   )]f6f|  cscl" 

(4-22)  d-q     =   {-  cos  cp     sin   (X-X   )   63 

+  cos  cp  sin  (X-X   )  [sin  (cp-cp  )  sin  (cp+cp  )  f5f3 


+  j-  tan  cp  sin  4  cp  f5f]  \   cscl" 


5a      2 
where     58  =  —  +  sin  cp  5  f 
a       Yo 

The  formula  (4-21)  is  quite  complicated,  but  can  be 
simplified  for  systems  not  of  world-wide  extent.  When  the  terms  f6f 
and  f5B  are  assumed  to  have  a  magnitude  of  1^6,000,000  or  less,  the 
terms  containing  the  squares  and  products  of  (cp-cp  )  and  (X-X   )   may  be 


29 


neglected  in  the  factors  by  which  the  terms  are  multiplied.  When 
these  angles  are  such  that  the  distance  from  the  origin  is  600  km. s 
the  squares  and  products  amount  to  an  order  of  magnitude  of 
1:600,000,000  (33). 

The  simplified  formulas  are:- 

(4-23)  d|  =  | [sin  (cp-cp  )  -  2  cos  cp  sin  cp  sin2  ^  (>  -X  )]°  P 

-  k   cos  cp  cos  -  (cp+cp  )  sin  -  (cp-cp  )5  f 

-  (2  +  4-  tan  cp  sin  4cp  )  sin  (cp  -cp  )ib  L'l   cscl" 

(4-24)  dt]  =  [-  cos  cp  sin  (X-X  )&  p 

+  j-  sin  cp  sin  4  cp  sin  (\-\  )  fof]  cscl" 

4.3     The  Use  of  Gravimetric  Deflections  of  the  Vertical. 

Probably  the  most  accurate  method  of  determining  tne 
equatorial  radius  from  deflections  of  the  vertical,  involves  a  final 
correction  by  application  of  tne  gravimetrically  determined 
deflections  of  the  vertical  ("absolute"  deflections). 

The  following  analysis  could  have  been  given  at  the 
beginning  of  this  chapter,  but  it  was  felt  better  to  place  it  here, 
thus  keeping  the  complete  picture  intact. 

Rice  (25)  p.  2,  states  that  "a  fundamental  geodetic  problem  is 
to  arrive  at  a  reference  ellipsoid  of  such  dimen;  ions  and  having  such 
orientation  that  the  deflections  of  the  vertical  are  as  small  as 
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possible  in  the  aggregate." 

An  ideal  reference  ellipsoid  would  be  centered  at  the  center 
of  mass  of  the  earth  and  its  polar  axis  would  coincide  with  the  mean 
axis  of  rotation  of  the  earth.  This  is  what  Veis  (32)  p.  l6,  calls 
the  Terrestrial  Ellipsoid.  A  further  requirement  which  would  be 
placed  on  an  ideal  reference  ellipsoid  would  be  that  its  surface  best 
fit  the  surface  of  the  geoid. 

The  origin  or  initial  point  of  a  geodetic  systerr  is 

generally  assigned  the  astronomic  latitude  and  longitude,  cp*  and  X', 

o     o 

observed  at  the  place.  The  geoid  distance  N*  is  generally  considered 
zero  and  no  errors  in  azimuth  are  assumed.  The  effect  of  making 
cp  =  cp',  and  \   =  X'at  the  initial  point  is  a  parallel  displacement  of 
the  reference  ellipsoid's  system  of  coordinates  with  respect  to  those 
of  the  terrestrial  ellipsoid  and  the  geoid.  At  another  point,  P.. ,  in 
the  same  geodetic  system,  which  is  displaced  from  the  geoid,  the 
normals  do  not  coincide,  the  result  not  only  of  the  displacement,  but 
of  the  wrong  computation  surface,  which  is  not  the  same  size  or  shape 
as  the  true  earth  spheroid.  The  resulting  meridian  angle  between  cp' 
and  cp  is  the  astro-geodetic  deflection  (cp*  -cp)  in  the  direction  of  the 
meridian.   (This  discussion  is  generally  in  terms  of  this  meridian 
component,  for  convenience,  though  similar  remarks  can  be  made 
concerning  the  prime  vertical  component,  as  well  as  the  total 
deflection). 

The  meridian  component  and  the  prime  vertical  component  of  the 
astro-geodetic  deflection  each  consists  of  three  angles  in  the  same 


31 


plane,  with  respect  to  the  normal  to  the  mean  or  smoothed  geoid.  If 
we  can  correct  this  deflection  in  such  a  manner  that  both  components 
are  reduced  to  &ero,  what  would  remain  would  be  the  normal  to  the 
mean  geoid .   In  doing  so  we  have  to  change  the  size,  shape,  and 
position  of  the  reference  or  computation  ellipsoid  so  that  our  point 
Pn  coincides  (as  nearly  as  possible)  with  the  corresponding  point  on 
the  mean  geoid,  and  so  the  new  ellipsoid  is  in  the  position  of  a 
terrestrial  ellipsoid.,  The  three  angles  ares- 

(1)  Correction  for  displacement,  (d|-y  <3\_) 

(2)  Correction  for  change  of  ellipsoid,  (d|  ,  dT  .J 

(3)  The  absolute  deflection  of  the  vertical,  or  the  angle 

between  the  true  vertical  and  the  normal  to  the 

mean  geoid,  or  if  preferred,  the  terrestrial 

ellipsoid  of  best  fit.  Since  this  is  obtained 

gravimetrically,  we  call  its  two  components  £ 

S 

and  ri  . 
g 

See  figures  3  and  k   for  pictorial  descriptions  of  this. 
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Fig.    k 
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It  should  be  noted  that  when  the  gravimetric  deflections  at 
the  initial  point  are  known,  two  of  the  unknowns  of  the  previous 
problem  are  now  known,  namely  5|  ,  and  §t]  .  These  are  equivalent  to 
the  following :- 

(4-25)     8  £  »  -  (S  -I)  =    U     -I) 
v  ''  bo     Vba   ^g'o   vsg   bayo 

6  ^o  =  "  K  '   Vo  =  {\  '   Vo 
where  (|a)Q  =  (q£  -  cpQ) 

^a}o  =  (Xo  "  \>>  C0S  ^o 

The  gravimetrically  computed  quantities  £  and  tj  ,  are  in- 

6      8 

dependent  of  the  size  of  the  reference  ellipsoid  (13)  p.  239;  even 
though  they  depend  on  gravity  anomalies  derived  from  the  gravity 
formula  employed. 

In  the  method  (described  in  section  4.1)  of  determining  the 

dimensions  of  the  earth  ellipsoid  by  means  of  correcting  the  astro- 

2 
geodetic  deflections  with  d|  and  d£  ,  and  subsequently  reducing  £| 

to  a  minimum,  we  have  one  effective  method  of  solving  the  problem. 

However,  we  can  go  one  step  farther,  as  indicated  in  the  preceding 

paragraphs,  and  make  the  final  correction  by  subtracting  the  value  | 

6 

from  the  astro-geodetic  deflections  which  have  already  been  corrected 
by  d£  and  d£  . 

Thus  our  problem  will  be  solved  by  a  least  squares  solution  in 
which 


3^ 


2 
£  (£  -  5  )  =  minimum 
6 

p 
2  (t)  -  T)  )  =  minimum 


In  this  way  we  are  much  closer  to  the  mean  geoid  at  the  start, 
The  result  should  be  a  terrestrial  ellipsoid  which  best  fits  the 
surface  of  the  geoid,  providing,  of  course,  that  we  have  sufficient 
points  from  a  large  enough  area  for  our  computation.  Theoretically 
the  problem  can  be  solved  when  astro -geodetic  and  gravimetric 
deflections  are  known  at  just  two  points,  one  at  each  end  of  a  long 
arc,  preferably  as  long  as  the  earth's  radius.  Naturally,  the  more 
points  at  which  this  data  is  known  the  better  the  solution.  Heiskanen 
has  suggested  this  method  in  a  number  of  publications;  see  par- 
ticularly (12),  and  (13)  p.  278-279- 


Chapter  5»  Flattening  from  Gravity  Anomalies. 

A  value  for  flattening  can  be  gotten  from  the  solution  by  the 
astro -gravimetric  method  described  in  section  4.3,  but  if  this  is  done 
the  gravity  anomalies  used  in  computation  of  the  gravimetric  deflect- 
ions are  no  longer  correct  since  these  anomalies  are  based  on  the 
adopted  gravity  formula,  which  in  turn  is  based  on  a  particular 
ellipsoid  with  its  own  flattening.  This  gives  a  clue  as  to  how  the 
flattening  may  be  corrected  frpm  the  anomalies  themselves.   See  (13), 

p.  279- 

According  to  the  theorem  of  Clairaut  (13)  p.  52,  the 
following  is  true  (extended  form):- 

(5-1)    p  =  i±rA  =  t m  - f  -  S  ■* 

E 

where    7    ,  7%   are  normal  gravity  at  pole,  equator,  resp. 

P  a>2 

m  =  (centrifugal  force  at  equator )/7 -  =  -~ 

E   7E 

f  =  flattening 

The  basic  form  of  the  gravity  formula  giving  normal  gravity  7 
for  any  latitude  cp  on  the  ellipsoid  is  (neglecting  the  longitude  term): 

(5-2)  7   =  7e  (1  +  p  sin2cp  -  €  sin2  2  cp) 

5      12 
where       e  =  fr  ntf  -  w  f 

The  problem  is  to  correct  the  parameters  7™,  (3,  and  e,  to  the 

gravity  formula  where  these  parameters  have  been  computed  for  a 

particular  reference  ellipsoid.  Thus  if  p  and  e  are  corrected,  a 
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correction  to  f  is  automatically  found.  The  correction  to  e  (based  on 
currently  used  ellipsoids)  is  always  so  small,  we  may  neglect  its 
computation. 

Using,  for  example,  the  International  Ellipsoid,  the  following 
values  are  used  as  a  basis :- 

a  =  6  378  388  m.  f  =  1/297  =  0.0033670 

7    =  978.0^9  gal  m  =  1/288.36  =  0.003^78 

0  =  0.0052884  €  =  0.0000059 

For  obtaining  a  corrected  value  of  7,  we  haves- 

7'  =  7  +   A7 

p 
where  A  7=  x'  +  y'  sin"  cp 

in  which  x'  =  unknown  correction  to  7  ,  in  mgals. 

JS 

y'  =  unknown  correction  to  |3,  in  mgals. 
Before  correction,  the  gravity  anomalies,  based  on  the 
ellipsoid  used  and  isostatically  reduced,  are:- 

Ag  =  gQ  -  7 
After  correction  tbey  are:- 

v  =  gQ  -  7'  =  S0  "  (7  +  A7) 
(gQ_7)  -  A7  =  Ag  -  A7 

Then,  we  can  write  error  equations  for  each  point  where 
gravity  observations  have  been  made t- 

x#  +  y  sin  cp  -  Ag  =  -  v 

If  we  have  enough  such  stations,  spread  out  over  a  large 
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enough  area,  we  can  solve  the  problem  effectively  by  the  method  of 

least  squares  in  such  a  manner  thats- 

2 
Z  v  =  minimum. 

An  example  of  this  type  of  solution  made  by  Heiskanen,  in- 
volved using  anomalies  from  over  1500  squares  of  1  by  1  dimensions 
(13)  P«  77  •  The  resulting  parameters  were:- 

7E  =  978o 0^51 
3  =  0.0053026 
e  =  0.0000059 

The  flattening  corresponding  to  the  corrected  £  was:- 

f  =  1/298.2. 


Chapter  6.  Gravimetric  Deflections  of  the  Vertical 

6.1  General. 

The  gravimetric  method,  of  which  the  subject  of  Chapter  5  is  a 
part,  also  produces  what  are  sometimes  referred  to  as  "absolute" 
deflections  of  the  vertical,  obtained  from  gravity  anomalies.  The 
absolute  or  gravimetric  deflection  of  the  vertical  at  a  point  is 
defined  as  the  angle  between  the  normal  to  the  geoid  and  the  normal  to 
the  earth  spheroid  at  the  point.   —  Veis  (32)  uses  the  term  "mean 
geoid"  to  describe  the  surface  of  the  earth  spheroid  in  the  area  of  a 
point,  meaning  a  smooth  surface  striking  a  "mean"  between  undulations 
of  the  geoid.   —  The  discussion  in  Chapter  k   (section  h.3),   describes 
the  use  of  these  deflections  in  determining  the  size  of  the  earth 
ellipsoid.  This  chapter  will  present  a  brief  outline  of  the  method  of 
determining  the  deflections  from  gravity  anomalies. 

To  begin  with,  the  gravity  anomalies  £g,  are  caused  by  dis- 
turbing mass  anomalies,  and  with  this  in  mind  and  with  the  aid  of 
fig.  5,  we  can  easily  see  how  local  absolute  deflections  occur. 


Fig.  5 
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The  mass  surplus  tends  to  "pull"  the  plumb  line  toward  it,  and 
conversely  the  mass  deficiency  tends  to  repel  it,  thus  causing  the 
deflection  i  =   cp'  -  cp.  The  sign  of  this  deflection  is  apparent  from 
this  expression  which  shows  that  |  is  always  referred  to  the 
astronomic  latitude,  cp! ,  of  the  point  in  question.  The  sign  conven- 
tion for  gravimetric  deflections  is  the  same  as  for  astro-geodetic 
deflections.  Referring  to  figure  5,   we  wee  that  when  cp  >  cp1 ,  £  = 
(cp'  -  cp)  <  0,  and  when  cp  <  cp1 ,  |  =  (cp'  -  cp)  >  0.   In  U.S.  practice  X 
is  considered  positive  westward^  and  in  European  practice  positive 
eastward.  In  either  case,  when  X  >  X' ,  t\  -   (X'  -  X)  cos  cp  <  0,  and 
when  X  <  X',  »  -  (X'  -  X)  cos  cp  >  0. 

From  a  theorem  developed  by  Stokes,  Vening  Meinesz  formulated 

expressions  for  absolute  deflections  of  the  vertical  in  1928.  Stokes' 

Formula  expressed  the  geoid  distances  N,  as  follows :- 

2rr    it 


(6-1)   N  =  2~~^U   A  J  f  (+)  sin  f  ^g  d  f 


in  which  \|r  =  angular  distance  of  a  circular  ring  of  width  d\Jr.     See 
fig.   6. 
Ag  =  average  gravity  anomaly  in  ring,  having  the  limits  \|/  and 
\|r  +  d\|r. 
R  =  Mean  radius  of  the  earth. 
G  =  Mean  gravity. 
f(+)  =  |  S(lr). 
(6-2)  S(\|r)  =  esc  s>  y  +  1  «  6  sin  ^  t  -  5  cos  \|r  - 


kO 


1       2  1 
-  3  cos  |ln  (sin  rf  +  sin  -^  t) 


which  is  the  original  Stokes  Function, 
sin  t  ch|r  dA  =  surface  area  element  for  unit  sphere, 

POLE 


P  —     Compu+a+ion 
Po/nt 

C-     =      Cen-f&r-       of 
Sphere, 

dS  =      ■Sw-ftkce.    Area. 
£  f  e  m  en  i~ 


Fig.  6 


Vening  Meinesz  got  his  formulas  by  differentiating  the  Stokes 
formula,  as  indicated  in  fig.  7 ,   in  the  form:- 

1 


(6-3) 


N  = 


2  7T  G  R 


f(\|f)  Ag  d  S 


GBOiD 


SPHtLISOlO 


r* 

<         dx 

• 

{ 

\ 

r 

^ — -~--£— 

dN 

Fig.    7 

i 

A/ OB.TH 


Upon  differentiating,  the  following  are  obtained:- 


in 


(6-k)     d|"  -  -  H  *  ^°  ^'g  /  f '  (♦)  cos  A  ^  d  S 

0  X   2  TT  G  IT  J 
s 

(6-5)    an"  -  -  H  -  -^^  f  f '  (♦)  sin  A  ^g  d  S 

y   2  f  G  R  j 
s 

2 
where     dS  =  R  sin  \|r  d  \|r  dA,  and 

cos  h 
(6-6)f '  (Mr)  =  ^M  =  -  °0S  |  jr  +  8  sin  *  -  6  cos  |* 
v       2  sin  ^ 

-  3(1  -  sin  ^)  esc  i|r  +  3  sin  t  In  [sin  ^|r  (1+  sin  ^)] 

Application  of  these  formulas  to  practical  use  requires 
further  manipulation.  They  are  used  in  different  form,  whether  one  is 
going  to  compute  the  effect  of  areas  beyond  t  ~  3  >   or  within  the  area 
circumscribed  by  a  circle  of  radius  \|r  ~  3  (3  is  an  arbitrary  dividing 
line,  which,  though  close  to  the  limit,  has  been  chosen  for  convenience 
(31)  p.  96).  These  different  forms  are  described  as  the  "Circle -Ring" 
Method,  for  use  up  to  \|r  «=  3  y   and  the  "Squares"  Method,  for  use 
beyond.  Still  another  form  is  used  for  the  small  inner  circle  close  to 
the  computation  point,  having  a  radius  r  ,  generally,  of  5  km.  or  less. 
Figure  8  shows  this  division  of  methods  pictorially.  For  further 
references,  see  (13)  p.  63 -70,  257-264,  (15). 

6.2  The  Squares  Method. 

This  method,  used  for  the  earth  represented  as  a  spherical 
surface,  and  for  the  area  beyond  a  radius  of  about  3  from  the  comput- 
ation point,  makes  use  of  mean  gravity  anomaly  values  for  each  1  by 
1  square  (spherical  trapezoid) between  \Jr  «=  3  and  \|r  ~  20  ,  and  for 
each  5  by  5  square  beyond  20  and  around  to  the  antipode  (31 )• 
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JZ    -     3°  x  S°     S&UAZES 


ig.  8 


These  should  be  either  free -air  or  isostatic  anomalies,  the  latter 
being  best.   If  Bouguer  anomalies  are  used  they  must  be  corrected  to 

obtain  free -air  anomalies.  To  adapt  formula  (6-k)   to  the  squares 

2 
method,  we  first  replace  dS  by  R  sin  \|r  d  \|r  d  A,  obtaining:- 


(6-t)      as"  = 


csc  1"    f 

2   TT   G      J 


f '  (\|r)  cos  A  Ag  sin  \|r  d  \|r  d  A 


Then  using  the  square  element  dq  on  the  unit  sphere,  R  =  1, 
(1*0,  we  obtain: - 

(6-8)   d|"  =  §2Lj£  /  f  MO  cos  A  Ag  d  q 

The  expression  for  &r\"   being  the  same  with  sin  A  replacing 
cos  A. 

o 

When  Z^  is  in  milligals,  and  when  G  is  taken  as  979*8  cm. /sec, 
the  term  before  the  integral  sign  will  be: 


^3 


c"  =  frT  =   0'.'0335 
One  procedure  is  to  perform  numerical  integration  summing  up 
the  effect  of  each  square  over  the  earth  as  indicated  in  the  following 
formulas  - 

(6-9)   d£"  =  c"Zf'0lOn  cos  A  £g  q 

o    o    o    o 

where  q  is  the  area  of  a  1  by  1  or  5  by  5  square  and  depends  upon 

its  mean  latitude.  The  subscript  q,  wher":  occurring,  signifies  mean 
values  of  the  quantity  for  the  square  in  question. 

The  squares  method  adapts  readily  to  solution  by  high  speed 
electronic  computing  machine,  but  manual  graphical  solution,  though 
more  time  consuming,  is  made  easier  by  use  of  maps  on  polar  stereo- 
graphic  projection.  With  this  projection  templates  can  be  easily 
prepared  for  use  on  the  large  area,  small  scale  maps  required,  since  a 
feature  of  the  projection  is  that  all  circles  on  the  sphere  project  as 
circles  on  the  map,  so  azimuth  lines  (which  are  great  circles)  and  the 
circular  rings  (small  circles)  can  be  easily  drawn.  Necessary 
formulas  for  preparing  such  templates  may  be  found  in  (13)  p.  259-260, 
and  in  (31)  p.  98-100. 

The  functions  f'OlO,  f'OlO  sin  ty,  and  /f'(\|r)  sin  f  d  t,  have 
been  computed  by  Sollins  at  close  intervals  for  different  angular 
distances  \|r,  from  the  computation  point.  These  tables  have  been 
published  by  Sollins  (29)  and  are  very  useful  in  the  solution  of  the 
problem. 


6» 3  The  Circle -Ring  Method. 

In  the  relatively  close  area  around  the  computation  point 
(within  \|r  ~  3  )>  where  a  detailed  gravity  anomaly  map  is  available, 

formula  (6-7)  is  used  in  the  following  form:- 

Ap  \|r+di|r 

(6-10)  V     f  f 

dr  =  ¥wir  J  cos  d  A  J  f'(t)  sin  *  %  d  * 

\|r+d\|; 
=  1*1  g"    (sin  A     -  sin  A^)  J    f'0lO  sin  «r  Ag     d  f 

For  d  T)",  substitute  (cos  A.  -  cos  Ap)  for  (sin  A  -  sin  A_). 

Formula  (6-10)  gives  the  effect  of  a  circle -ring  compartment  bounded 

by  ty  and  \Jr  +  di|r  and  by  A.  and  A^,  and  having  a  mean  anomaly  ££,>   on 

the  total  component  £  ". 

r  g 
Values  of  J~  "''(t)  sin  ty  d  \|r  are  tabulated  in  Sollins  Tables. 
^=10m  \|r+dilr 

Thus  for  the  value  of  a  particular  /f'  (t)  sin  \|r  d  \|r,   the  tables 

* 

would  be  used  as  indicated  in  the  following  expression: - 

\|r+d\|r  \|r 

/  f  (\|r)  sin  t  d  t  -  /  f'MO  sin  ♦  d  t. 

t=10m  \|r=10m 

Figure  9  gives  pictorial  explanation  of  the  components  of  the 
method  and  corresponding  templates. 

Any  system  of  azimuth  sector  and  ring  zone  intervals  can  be 
used  in  this  method,  but  a  system  which  yields  equal  effect  for  each 
compartment  would  obviate  much  unnecessary  work.  Kasansky  developed 
a  system  in  which  the  sine  differences  are  l/8  (32  compartments),  and 


A/O&TH 


^ 


R.TME.MT 


P  =  computation  point 
\(f  =  inner  radius  of  compartment 
\|r  +  d\|r  =  outer  radius  of  compartment 
A. ,  Ap  =  bounding  azimuth  lines 
A  g.  =  mean  gravity  anomaly  in  mgals,  read  from 
is o -anomaly  lines  on  map 
Fig.  9.  Components  of  Template  for  Circle - 
Ring  Method 
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the  ring  radii,  beginning  with  r  =  5.0  km.,  follow  the  geometrical 
progression  r  =  1.270  r  .  Farther  from  the  station,  the  ratio  is 
slightly  smaller  than  I. 270  because  of  the  earth* s  curvature.  Tables 
of  these  Kasansky  Circle -Rings  may  be  found  in  (13)  p.  267,  and  in 
(16)  p„  17.  The  effect  of  each  compartment  is  O'.'OOl  radial  deflection 
effect  per  1  mgal  Ag. 

D«  A.  Rice,  of  U.S.C.&  G.S.,  published  a  system  of  Circle - 
Rings  (25)5  in  which  each  ring  was  divided  into  36  equal  azimuth 
sectors  of  10  aperture.  Here  also,  the  radial  effect  for  each  com- 
partment for  1  mgal  Ag  is  O'.'OOl.  Rice's  method  for  computation  of  the 
radii  involved  using  the  formula  for  the  deflection  effect  along  any 
radial  lines- 


csc  1"     /      , 
2  ?r  G     J 


(6-11)  d0"    ::    ^    ^       J      f»0|r)    Sin    t    d   *    Ag    d    A 

s 


so  d£"  =  d  9"   cos  A 

dn"  =  d  0"   sin  A,   and 

(6-12)        d0"   =  |2S^1     Ag  ijl  M 

where  I.   is  obtained  from  Sollins'  Tables  (29). 
*o 

Then  for  AA  =  10°  =  0. 17^52  radians,  and  for  G  =  981*000  mgals, 

formula  (6-11)  becomes?- 

*1 

d0"  =  0.0058^05^8^  I  x   Ag. 

Vo 
For  d0"  =  O'.'OOl  and  Ag  =  1  mgal, 

AX     =    X  ^  -   0.17121680 
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*2   +1 
AI2  =  I.   -  I   =  0.17121680,  etc. 

o     o 
Then  selecting  interpolated  values  of  \Jr  from  Sollins1  table 
giving  AI,  the  circle -ring  table  was  prepared.  Table  I  is  Rice's 
table  reproduced. 

Table  I 

Rice  Template  Radii  for  Gravimetric  Deflection  of  the  Vertical 

For  radial  deflection  effect  of  O'.'OOl  'and  mean  com- 

s  o 

partment  anomaly  of  1  mgal;  angular  aperture  10  . 

(n  =  zone  number;  r  =  inner  radius  in  km. ) 
n       rn       rn       rn       r 


1 

0.119 

14 

1.099 

27 

10.15 

39 

77-97 

2 

O.lfcl 

15 

1.304 

28 

12.05 

40 

92.22 

3 

O.167 

16 

1.5V7 

29 

14.29 

41 

109.0 

4 

0.198 

17 

I.836 

30 

16.94 

42 

128.7 

5 

0.235 

18 

2.179 

31 

20.09 

43 

151.9 

6 

0.279 

19 

2.586 

32 

23.83 

44 

179.1 

7 

0.331 

20 

3.068 

33 

28.25 

45 

210.9 

8 

0.393 

21 

3.641 

34 

33.48 

46 

248.0 

9 

0.467 

22 

4.320 

35 

39.67 

47 

291.2 

10 

0.554 

23 

5.125 

36 

47.00 

48 

341.2 

11 

0.657 

24 

6.081 

37 

55.66 

^9 

399.0 

12 

O.780 

25 

7.216 

38 

65.90 

50 

465.5 

13 

O.926 

26 

8.560 

With  the  use  of  circle  rings  such  as  Rice's,  the  problem  of 
computing  the  gravimetric  deflections  becomes  simple,  though  laborious. 
All  that  is  necessary  is  to  estimate  the  mean  anomaly  from  the  contours 
in  each  compartment,  make  summation  of  Ag  for  each  azimuth  sector,  mul- 
tiply by  0.001  cos  A  or  0.001  sin  A  ,  and  the  effect  for  that  sector  is 

m  m 

solved.  Then  the  effects  of  all  36  sectors  are  added  together  for  total 
effect.  This  may  be  further  simplified  when  one  takes  note  of 


i+8 


the  fact  that  mean  azimuths  in  ltfO  opposite  sectors  will  have  the 
same  cosine  or  sine,  but  with  opposite  sign.  One  must  be  careful  with 
signs,  since  both  the  azimuth  of  the  sector  and  the  anomaly  gradient 
contribute  to  the  sign.  Additional  information  on  this  aspect  of  the 
problem  will  be  given  in  the  explanation  accompanying  the  computations 
made  in  Chapter  7« 

In  the  event  gravity  material  is  poor,  which  is  too  often  the 
case  (and  one  has  to  be  a  good  judge  concerning  when  poor  material  is 
to  be  considered  good),  Heiskanen  (13)  p.  268,  recommends  combining 
two  zones  and  two  sectors  into  one  four-compartment  compartment,  the 
effect  of  which  then  would  be  O'.'OOU  per  mgal. 

6.k     The  Inner  Circle. 

The  circle -ring  method,  although  it  can  be  applied  beyond 
ty  =  3  y   should  not  be  used  very  close  to  the  computation  point  since 
the  term  esc  -  ty  in  the  value  of  f'(\|r)  sin  \|r,  in  formula  (6-10) 
approaches  infinity.  Rice  (25),  suggests  a  lower  limit  of  100  meters. 
However  this  is  usually  too  small  because  of  the  spars ity  of  gravity 
anomaly  values,  except  in  isolated  instances,  within  any  given  100 
meter  circle. 

The  inner  circle,  being  small,  is  considered  to  have  a 
uniform  anomaly  gradient,  so  the  anomalies  read  from  the  cardinal 
points  are  used  in  the  following  expressions,  the  derivation  of  which 
may  be  found  in  (31)  p.  93,  100-101:- 

(6.13)   ar  .  £|£^  L&l  x  r  .  0V105  rQ  !i£L* 


h9 


2G     &y     o      ^  o   8  y 

where       r  -  radius  of  inner  circle  in  km. 
o 

5  (^€)x  =  ^€s  -  £gN  in  mgals  (N.  hemis. ) 

6  (£g)   =  Z&  -  £g       in  mgals  (European  practice) 

ox      =  2  r  ,        5y    =  2  r 
o  o 

Substituting  these  where  appropriate,  formulas  (6-13)  become: - 

(6-14)   a$"  =  0V0525  (^gg  *-  AgN) 
at)"  =  OV0525  (£xr   -  AgJ 

W       £j 

See  figure  10  for  further  clarification  of  this  "one- 

gradient"  methoa  of  computing  the  effect  of  the  inner  circle. 

x 


A?J =+' \^6e — *-y 


Fig.  10.  Inner  Circle  -  One -Gradient  Methoa 

Rice  usea  a  three -gradient  methoa  rather  than  one -gradient  as 
above,  in  computations  for  l6  points  in  the  U.S.  (26).  Taking 
anomalies  at  NE,  SE,  SW,  ana  NW  points,  as  well  as  at  the  carainal 
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points  of  the  circle,  three  N-S  gradients  and  three  E-W  gradients  were 
used,  with  weights  of  one -half  given  to  the  outer  gradients  and  a 
weight  of  one  to  the  center  gradient.  Figure  11  illustrates  this. 
This  method  should  be  used  only  when  detailed  gravity  information  is 
available  inside  the  inner  circle. 


*3w 


hx  rs-A/) 

<5x  (Sw-s/vJ) 
(SE-NE) 

Stj      (Sw-SE) 
(A/W-NE) 


=  2  r. 


>  —     1.4-14-   /"„ 


Fig.  11.   Inner  Circle  -  Three -Gradient  Method 
4g„  -  ^ 


d|"(S-N)   =  0'.'105 


'N 


8  x 


Weight  =  1 


d|"(SE-NE)  =  0V105  —^ —  r     Weight  =  ± 

O  X        O  l_ 

d|*'(SW-NW)  =  0V105   SW6  x  m   rQ    Weight  =  | 

£  Weights  =  2 

Combining  these,  using  weights  shown  and  values  of  ox  given 
in  fig.  11,  we  obtain  the  following  expression,  given  in  (31)  p.  101: 


(6-15)   d|"  =  0.0262  (Ags-  4sN)  +  0.0186  (^-z^+z^-Ag^) 
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This  is  for  the  northern  hemisphere.  For  the  southern 
hemisphere,  just  reverse  S  and  N  where  occurring.  The  other  component 
of  the  deflection  is  similars- 

(6-16)  m"  =  0.0262  (ZV^gjg)  +  °-0186   ^sw^sE^"  ^m] 

This  is  for  European  practice,  in  which  K   and  t]  are  considered 
positive  eastward.  Where  the  practice  is  opposite,  W  and  E  should  be 
reversed  where  occurring. 

As  may  be  deduced  from  the  opening  paragraph,  the  smaller  the 
inner  circle,  down  to  certain  limits,  the  more  concentrated  the  gravity 
information  must  be.  Vening  Meinesz  (13)  p.  70,  states  that  the  error 
incurred  using  a  circle  with  r  =  30  km.   is  less  than  1$.  A  good 
round  number  seems  to  be  5  km. ,  and  for  good  results  there  should  be 
from  four  to  eight  points  of  known  Z^  within  the  circle  (13)  p.  267. 


Chapter  '(.      Computations 

7.1  Dimensions  from  the  98th  Meridian  in  the  United  States. 

(From  Astro-Geodetic  Deflections). 

Using  data  from  U.S.C.  &  G.S.  Special  Publication  No.  229,  (7), 
consisting  of  astronomic  and  geodetic  latitudes  and  meridian  com- 
ponents of  the  astro-geodetic  deflections  of  the  vertical  for  30 
points  between  97  W  and  99  W  longitude,  extending  from  North  Dakota 
to  Texas,  values  for  a  and  f  will  be  determined.  The  astro-geodetic 
deflections  as  published  in  the  publication  named  have  been  reduced 
according  to  the  Hayford  isostatic  assumption  for  a  depth  of  com- 
pensation of  113  •  7  ton. 

2 
The  problem  to  be  solved  is  £  %     =  minimum,  according  to  the 

theory  of  least  squares.  The  Vening  Meinesz  displacement  and  change  of 
ellipsoid  formulas  will  be  used.  To  adapt  them  to  practical  use,  terms 
containing  the  same  unknowns  must  be  collected  and  combined. 

Prior  to  solution,  a  test  problem  was  solved  to  determine  what 
terms  may  be  disregarded  and  whether  or  not  the  long  or  short  formulas 
should  be  used.  The  maximum  difference  in  longitude  is  about  1  ,  and 
terms  containing  longitude  were  determined  to  make  negligible  con- 
tribution to  the  residual  deflection.  Also,  for  the  arcs  of  limited 
length  (A  cp  ~  13^3*  in  this  case)  encountered  in  this  problem,  the 
difference  between  -residuals  computed  by  long  and  short  formulas, 
turns  out  to  be  negligible.  Terms  containing  8N  will  also  be 
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neglected  to  simplify  the  problem. 

The  Vening  Meinesz  formulas,  given  as  (4-19)  and  (4-23)  in 
section  4.24,  are  here  repeated  with  terms  omitted  as  noted  above. 

Displacements-       d£  =  cos  (cp  -  cp  )  6| 

Change  of  Ellipsoids-  dg  =  P-.&P  -  p»6f  -  p~f  5  f 

Ba      2 
where  6P  =  —  +  sin  cp  6  f 

3.  O 

Substituting  the  foregoing  expression  for  5£  in  d£  and  collecting 

terms,  we  obtains - 

Sa        2 
<U  =  P-,  —  +  P-,  sin  cp  5f  -  p0&f  -  p„f&f 
*c   *±  a   ^1     ^o     r2  3 

pl  2 

=  —  5a  +  [Pl  sin  cpo  -  p2  -  p3f  ]  5f 

in  which  p.  =  sin  (cp  -  cp  ) 
*l      VYn  ^oJ 

p0  =  4  cos  cp  cos  ^(cp  +cp  )  sin  ^-(cp  -  cp  ) 
*2  Yn    2VYn  Yo'     2  n  Yo' 

P3  =  (2  +  ^  tancpQsin  1kpq)  sin  (cpQ  -  cpQ) 

The  observation  equations  will  take  the  form  given  below,  with  sub- 
script 0  signifying  the  initial  point,  Meades  Ranch,  Kansas,  and  sub- 
script n  signifying  any  other  point  in  the  arc. 

|  =  d|_  +  d£  +  (%%    -  cp  ) 
bn    bD    bc   N  n   Yny 

For  convenience  we  transform  the  unknowns  and  the  coefficients 

as  indicated  heres- 

5  a 


in   =  (h)S|0  -  (J)  1000  +  (k)  10,000  Sf  +  (1) 
=  hx  +  jy  +  kz  +  1 


^ 


where       x  =  5| 


and 


y  - 

z  = 

h  - 

0  = 

k  = 


o 
5a 


1000 
10,000  8f 


COS     (CD       -    CD    ) 

VYn       Yo/ 


sin   (cp     -  cp  ) 

VYn       Yo' 


>     unknowns . 


1000  esc  1" 


[sin   (cpn-<P0)[   sin2  cpQ  -   (2+|  tan  cpQ  sin  k  cpQ)f] 

I  1    /       .       \       •      1    /  \1   esc   1" 

-  k  cos^cp     cos  -  (cp  +cp   )  sin  =■  (cp  -cp  )  ^^^ 

^n  2   VYn  Yo'  2    v^n  Yo'J  10,000 


JJ 


Table  II  -  Basic 

Data  -  98th 

1  Meridian 

Point 

Latitude 

Astro -Geodetic 

n 

Location 

^r, 

Deflection 

n 

(cp'  -  cp) 

i+96 

N.Dak. 

l46°02'17"N. 

IV 17 

571 

S.Dak. 

45  28  16 

1.06 

566 

11 

44  53  53 

0.78 

569 

11 

44  00  52 

-2.93 

570 

n 

43  43  24 

-1.43 

567 

M 

>3  18  48 

O.96 

412 

Nebr. 

42  25  29 

3.88 

407 

11 

41  35  44 

-1.29 

413 

M 

4o  46  35 

-2.93 

235 

Kans . 

39  46  35 

-2.03 

Meades  R» 

11 

39  13  27 

0.00 

238 

11 

38  43  47 

-1.29 

2U5 

11 

37  51  38 

-2.51 

2i4-0 

11 

37  14  28 

-0.40 

522 

Okla. 

36  4o  08 

0.20 

524 

11 

35  56  49 

1.88 

516 

11 

35  16  27 

1.71 

521 

11 

34  38  53 

1.98 

602 

Texas 

33  37  21 

-1.28 

620 

M 

32  48  10 

-0.91 

639 

11 

32  15  45 

6.78 

605 

M 

31  39  26 

2.20 

594 

11 

30  59  k6 

7.65 

592 

11 

30  16  21 

-1.11 

635 

•1 

29  42  49 

-4.05 

652 

11 

29  22  16 

2.14 

627 

11 

28  52  38 

-0.59 

634 

ti 

28  18  06 

-0.68 

588 

11 

27  44  35 

1.44 

630 

11 

27  4o  16 

4.91 

Clarke's 

1866  Ellipsoid:. 

a  =  6, 

378,206 

m. 

f  =  1/294.98 

=  0. 

00339006 

Initial 

Points-   Meades 

Ranch 

q>0  =  39°13'26V686  n. 
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Table  III  -  Observation  Equations  -  98th  Meridian 


Point 

h 

J 

k 

n 

(x) 

(y) 

00 

=  (cp'    -cp  ) 

496 

0.993 

3.837 

-1.544 

IV 17 

571 

.99^ 

3.519 

-1.450 

1.06 

566 

.995 

3.197 

-1.348 

0.78 

569 

•  997 

2.701 

-1.178 

-2.93 

570 

•  997 

2.535 

-1.118 

-1.43 

567 

•  997 

2.306 

-1.033 

O.96 

412 

.998  • 

1.806 

-0.835 

3.88 

4o7 

.999 

1.338 

-  .637 

-1.29 

H3 

1.000 

0.876 

-  .429 

-2.93 

235 

1.000 

.312 

-  .158 

-2.03 

Meade s  Ranch 

1.000 

-- 

-- 

0.00 

238 

1.000 

-  .279 

.146 

-1.29 

245 

1.000 

-  .770 

.414 

-2.51 

240 

.999 

-1.119 

.613 

-0.40 

522 

.999 

-1.442 

.803 

0.20 

52  4 

.998 

-1.849 

1.051 

1.88 

516 

.998 

-2.228 

I.290 

1.71 

521 

.997 

-2.580 

1.520 

I.98 

602 

.995 

-3.157 

1.911 

-1.28 

620 

.994 

-3.617 

2.367 

-0.91 

639 

.993 

-3.920 

2.457 

6.78 

605 

.991 

-4.258 

2.710 

2.20 

594 

.990 

-4.628 

2.993 

7.65 

592 

.988 

-5.032 

3.310 

-1.11 

635 

.986 

-5.353 

3.567 

-4.05 

652 

.985 

-5.534 

3.716 

2.14 

627 

.984 

-5.808 

3.944 

-0.59 

634 

.982 

-6.128 

4.214 

-0.68 

588 

.980 

-6.437 

4.481 

1.44 

630 

.980 

-6.477 

4.516 

4.91 
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Table  IV  -  Least  Squares  Solution  -  98th  Meridian 


Forward  Solution; - 

x        y        z        1        L 

x  +29.62     -i'7.43         +35-76         +15  06  +33-03 

-1        +1.601263  -1.207292  -0.50911^  -1.115125 

y  +398-56        -250.22        -88.55  +12.36 

-  75-95      +  57.26      +24.15  +52.89 

(1)  +322.61      -192.96       -64.40  +65.25 

-1  +0.598122  +0.199622  -0.202257 

z  +160.48       +58,88  +  4.90 

-  43.17       -18.21  -39.88 
-115.41       -38.52  +39.03 

(    )  +     1.90       +  2.15  +  4.05 

-  1  -1.131579  -2.131579 

1  +223.62  +209.03 

-  7.68  -  16.82 

-  12.86  +  13.03 

-  2.43  -  U.58 

[||]  +200.65  +200.66 

Back  Solution; - 

z       Check       y      Check      x      Check 

1  -1.131579  -2.131579  +0.199622  -0.202257  -0.509115  -1.115125 
z  -O.676822  -1.274944  +I.366 146  +2.573438 
y  -O.764I32  -2.365417 

2  -1.131579  -2.131579  -0.U77200  -1.47720.1  +0.092899  -0.907104 
Z(check-U)         -1.000000          -1.000001  -1.000003 


Standard  Error  of  single  observation;- 


m  =  *fW-  +-f^^-  t-rnzm-wi* 


58 


Table  V     -     Computation  of  Correlate  Numbers   -  98th  Meridian 

Kx)  l(y)  l(z) 

(x)  +  1 

(y)  0 

+1.6013 

+1.6013  +  1 

(z) 


0 

0 

-1.2073 

,0 

+0.9578 

+0.5981 

-0.2495  +0.5981  +  1 

i-9o  o!zz  =  1  1.90  q     =  0.5981     1.90  a.  =  -0.2495 

Qzz  =  O.5263  Q^  =  0.3148  Qxz  =  -0.1313 

322.61  Q^  -  192.96  Q^  =  1 

322.61  o=i  +  60.7^38  =  61.7438 

0^=0.1914 

322.61   ft       -   192.96    Qa;x  =    1.6013 
322.61  G^x  =  1.6013   -  25.3356  =   -  23.7343 
0^  =  -  0.07357 

29.62  Q^  -  47.43  Q^  +  35.76  Qxz  =  1 
29.62  Qjqj  =  1  -  3.^89^  +  4.6953  =  2.2059 
Q^  =  0.07447 


Standard  Errors  of  Unknowns  after  Adjustment :- 

m    =  |W   0^       =     2.726  n/  0.67447  =  2.726   (+  0.2729)  =  +  0.7439 

«       =     lW     0^  =        2.726   yj     0.1914       =     2.726      (+    0.4375)     =     +    1.1926 


m    = 
z 


=     llsj     Q^  =        2.726   s]     0.5263       =    2.726     (+    0.2294)    =    +    0.6253 
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5a 


Results  of  least  squares  solution: - 
x  =  +  O.092899   =  55 

y  =  -  0.477200    1000 

z  =  -  1.131579   =  10,000  5f 
8 1  -  +  0709  +  0"lk 
5a  =  1000  (-  0.4772)  =  -  ^77.2  m. 

5f  "  '  io*ooo79  -  -  o-oooii3i58 

a  =  6',  378, 206  m. 
5a  =      -477 
New  a  =  6,377,729  m.  +  1190  m. 

f  =  0.00339006o 

5f  =  -0.000113158 
New  f  =  0.003276902  +  0.00006253 
l/f  =  304.9  +5.4 

The  solution,  as  can  be  readily  seen,  does  not  conform  to  what 
is  known  to  be  the  truth,  and  the  standard  errors  resulting  show  that 
the  precision  of  the  solution  is  "somewhat"  rough.  This  of  course  can 
be  explained  by  noting  that  the  length  of  the  arc  is  quite  short,  and 
that  the  number  of  points  used  is  small.  Better  results  certainly 
could  have  been  obtained  if  many  points  having  greater  longitude 
differences  had  been  used,  thus  affording  the  benefit  of  the  prime 
vertical  components  of  the  deflections.  However,  the  method  used  is 
what  is  intended  to  be  demonstrated  here. 
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It  could  have  been  expected  that  the  results  obtained  from  the 
data  used  would  be  as  they  are.  Differences  in  latitude  between  the 
initial  point  and  other  points  in  the  arc  apparently  differ  systemat- 
ically.  It  is  seen  as  a  rough  check  of  the  standard  error  of  the 
value  of  a  obtained  in  the  least  squares  solution,  that  if  we  take  the 
average  value  of  (cp*  -qp)  for  the  five  northernmost  stations  and  of  the 
five  southernmost  stations,  the  resulting  difference  between  the  two 
is  -  8'.'57<>  For  the  arc  as  short  as  the  one  used  this  is  high. 

Taking  the  mean  latitudes  of  each  of  the  two  groups  of  five 
points,  we  get  Acp  =  16  26'.  Taking  the  length  of  A  (cp'-cp)  at  31 
meters  per  second  of  arc,  we  have  -  8'.'5T  @  31  m.  =  -266  meters.  Using 
one  radian  as  57.3*  we  can  determine  the  proportional  difference  for 
an  arc  of  the  length  of  the  earth's  radius :- 

Sl|  (  -  266)  =  -  930  meters. 

This  kind  of  rough  calculation  could  tell  us  in  advance 
approximately  what  results  can  be  expected. 
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To 2  Dimensions  from  Sixteen  Stations  in  South  Central  United  States 

by  the  Astro-Gravimetric  Method. 

For  an  illustration  of  the  astro-gravimetric  method  as 
described  in  section  *f.3>  the  published  data  on  l6  points  in  South 
Central  U.S.  (including  Kansas,  Missouri,  Oklahoma,  Arkansas,  Texas, 
and  Louisiana)  as  computed  by  Rice  (26),  was  used. 

In  the  referenced  publication,  all  geodetic  data  were  con- 
verted from  Clar^c',r-  H866  Ellipsoid  to  the  International  Ellipsoid.   In 
addition,  as  a  result  of  the  investigations,  Rice  determined  that 
several  of  the  astronomic  observations  which  had  been  published  in  (7) 
were  wrong,  and  were  consequently  done  over.  Moreover,  all  latitude 
and  longitude  determinations  were  reduced  to  the  FK~  system  of 
cataloging  mean  places  of  stars.  The  gravity  anomalies  used  were  free- 
air. 

All  the  data  needri  for  solution  of  the  present  problem  was 
given  in  Tables  VT  and  VTI  of  ref.  (26),  namely  the  latitudes  and 
longitudes  of  each  station,  as  well  as  the  astro -geodetic  and  the 
gravimetric  deflections  for  each.  Here  again,  as  in  the  problem  of 
section  7-lj  the  extent  of  the  area  covered,  and  the  number  of  stations 
involved  are  inadequate  for  a  meaningful  solution  to  the  problem,  but 
for  illustrative  purposes,  the  results  are  gratifying. 

Meade's  Ranch,  Kansas,  was  again  chosen  as  the  initial  point. 
The  Vening  Meinesz  abbreviated  formulas  were  determined  to  be  adequate 
and  were  used  for  both  £   and  t).  Since  Meade's  Ranch  was  one  of  the  16 
points  given,  the  values  of  b£     and  br\     are  known,  thus  reducing  the 
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number  of  unknowns  to  two  (5N  again  being  neglected). 

In  adapting  the  formulas  (14—19),  (^-23),  (^-20),  and  (k-2k)  to 
this  problem,  the  terms  were  grouped  in  the  logical  manner  for  ease  in 
computation  of  the  coefficients  as  follow: - 

|  =  d|D  +  d|q  +  (  cp'  -  cp)  -  £g 

6  =  gl5^Q  +  h15lQ  +  ^  ^jg  +  ^  (10,000)  6f  +  (  cp'  -  cp  )  -  £g 

gx  =  +  cos  (cp  -  cpQ) 

h.  =  +  sin  cp  sin  (X-X  ) 

J,  =  +  [sin  (cp  -  cp  )  -  2  cos  cp  sin  cp  sin2  ~  (X  -  \  )]  —  esc  1" 


kl-  + 


2  1  2 

[sin  (cp  -  cp  )  -  2  cos  cp  sin  cp  sin  ^  (X-X  )]  sin  cp 

-  k  cos  cp  cos  g-  (cp  +  cpQ)  sin  g-  (cp  -  cpQ) 

-  (2  +  |  tan  cpQ  sin  k  cpQ)  sin  (cp  -  <P0)f]g£-£ 


T)  =   dV  +  dTi      +    L  V    -  X)    COS   cp   -T] 

D  c  g 

ti  =  g28  £o  +  h15T)o  +  j2  —  +  k2   (10,000)   8f  +   (X«    -  X)   cos  cp  -  % 

g2  =   -  sin  cpQ  sin  (X  -  Xq) 
h2  =  +  1 

j2  =   -  cos  cpQ  sin   (X-  Xq)  ~2  Csc  1" 

k«  =   I    -  cos  cp     sin  (X  -  X   )  sin     cp 
2  o  x  o'  MD 


+  J  sin  cpQ  sin  k  cpQ  sin  (X-  \Q)  f]  ||£ 


csc  1" 


000 
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Then  the  observation  equations  "becomes- 

1  =  h  loo  +  ki  (10>00°)  6f  +  \ 

where  1±  =  g;L  8£q  +  h^  tjq  +   (|&  -  |g) 
ia  =  <P'   -  9 
6  iQ  =   -   (ia  -   ^g)0  =    (Ig  -  ia)Q 


n  =  j2  ig§  +  k2  (10,000)  8f  +  i2 

where  12  =  g25|Q  +  hg5  tjq  +   (t^  -  tj    ) 

na  =  (v  -  x)  cos  cp 
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Table  VI  -  Basic  Data 

-  Rice' 

s  16 

Points 

Point 

n 

Name 

? 

i 

\ 

1 

Meade s  R. 

39° 

13' 

27"  N. 

98° 

32 , 

31"  w 

2 

Arcadia 

37 

30 

08 

9k 

35 

kl 

3 

Dirks 

37 

25 

40 

99 

16 

21 

4 

Twin 

36 

06 

57 

96 

k6 

40 

5 

Little  Rock 

34 

44 

51 

92 

16 

25 

6 

Bogue 

33 

43 

05 

90 

53 

17 

7 

Ecore 

33 

36 

20 

92 

54 

25 

8 

Polk 

33 

23 

05 

99 

05 

46 

9 

Bartley 

4  33 

12 

08 

101 

ho 

44 

10 

Roby 

32 

kk 

36 

100 

22 

46 

11 

Lacasa 

32 

39 

05 

98 

hi 

30 

12 

Burns 

32 

35 

28 

93 

14 

20 

13 

Sears 

32 

33 

31 

100 

02 

16 

14 

Brooks 

32 

25 

09 

92 

07 

01 

15 

Bynum 

32 

19 

16 

101 

00 

00 

16 

Legion 

32 

16 

19 

90 

06 

07 

International  Ellipsoid:-   a  =  6,378,308  m. 

f  =  1/297 

=  0.003367003367 

Initial  Point:-     Meades  Ranch 

q>    =  39°  13'   26V686  N. 
X     =  98°  32"    30'.' 506  W. 


bo       Vbg       *a'o 


-  Ua  -  6g)0  =  -(-0V90)  =  +  0V90 


6  1      =    (T)      -  ti    )      = 
'o       v  g         a'o 


-   (\  -  Tlg)0  =   -(-0V30)  =  +  0V30 
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Table  VI     (cont'd) 


(1) 

(2) 

(3) 

(4) 

>oint 

gl 

5    *o 

hl 

6    \ 

5a-*g 

h 

(0*.' 

90  gx) 

(OV 

30)  h1 

11 

(1)  +  (2)  +  (3) 

1 

+ 

0.90 

0 

- 

0.90 

C 

2 

+ 

.90 

- 

0.01 

- 

.16 

+  O.'ii 

3 

+ 

•  90 

+ 

.00 

- 

.98 

-      .C*3 

4 

+ 

.90 

- 

.01 

- 

.73 

+      .16 

5 

+ 

.90 

- 

.02 

+ 

.66 

+  1.54 

6 

+ 

.90 

- 

.02 

+ 

.56 

+  1.44 

7 

+ 

.90 

- 

,.02 

+ 

.07 

+     .95 

8 

+ 

•  90 

+ 

.00 

- 

.75 

+     .15 

9 

+ 

.90 

+ 

.01 

- 

.89 

+      .02 

10 

+ 

.89 

+ 

.01 

- 

.36 

+      .54 

11 

+ 

.89 

+ 

.00 

- 

.25 

+     .64 

12 

+ 

.89 

- 

.01 

+ 

.56 

+  1.44 

13 

+ 

.89 

+ 

.00 

- 

.50 

+     .39 

14 

+ 

.89 

- 

.02 

- 

.06 

+     .81 

15 

+ 

.89 

+ 

.01 

- 

.45 

+     .45 

16 

+ 

.89 

- 

.02 

+ 

.44 

+  1.31 

taint 

g2 

5   *o 

h2 

5   \ 

^a 

im\ 

h 

(ov 

90  gg) 

(0'.'30  h2) 

11 

it 

(1)  +  (2)  +  (3 

1 

0 

+ 

0.30 

m 

0.30 

0 

2 

+ 

0.04 

+ 

.30 

- 

1.24 

-  0.90 

3 

- 

.01 

+ 

.30 

- 

.34 

-   .05 

k 

+ 

.02 

+ 

.30 

- 

.61 

-   .29 

5 

+ 

.06 

+ 

.30 

0.00 

+    .36 

6 

+ 

.08 

+ 

.30 

- 

.05 

+   .33 

7 

+ 

.06 

+ 

.30 

+ 

.76 

+  1.12 

8 

- 

.01 

+ 

.30 

- 

.02 

+      .27 

9 

- 

.03 

+ 

.30 

+ 

2.26 

+  2.53 

10 

- 

.02 

+ 

.30 

+ 

•  74 

+  1.02 

n 

- 

.00 

+ 

.30 

+ 

.85 

+  1.15 

12 

+ 

.05 

+ 

.30 

+ 

.25 

+      .60 

13 

- 

.01 

+ 

.30 

+ 

1.14 

+  1.43 

14 

+ 

.06 

+ 

.30 

+ 

.47 

+      .83 

15 

- 

.02 

+ 

.30 

+ 

1.88 

+  2.16 

16 

+ 

.08 

+ 

.30 

+ 

.27 

+     .65 
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Table  VII     -     Observation  Eauations      -     Rice's  l6  Points 


Point 


1 

0 

0 

0 

2 

- 

0.093 

+ 

0.672 

+     .73 

3 

- 

.102 

+ 

.773 

-     .08 

k 

- 

.176 

+ 

I.366 

+     .16 

5 

- 

o26l 

+ 

1.993 

+  1.54 

6 

- 

.323 

+ 

2.484 

+  1.44 

7 

- 

.323 

+ 

2.554 

+     .95 

8 

- 

.329 

+ 

2.681 

+     .15 

9 

- 

.341 

+ 

2.767 

+      .02 

10 

- 

.366 

+ 

3.001 

+     .54 

11 

- 

.370 

+ 

3.050 

+     .64 

12 

- 

.379 

+ 

3.065 

+  1.44 

13 

- 

0^76 

+ 

3.096 

+     .39 

14 

- 

.392 

+ 

3.1^ 

+     .81 

15 

- 

.390 

+ 

3.215 

+     .45 

16 

- 

.4o6 

+ 

3.205 

+  1.31 

1 

0 

0 

0 

2 

- 

0.172 

+ 

0.440 

-  0.90 

3 

+ 

.032 

- 

.081 

-     .05 

4 

- 

.077 

+ 

.197 

-     .29 

5 

- 

.275 

+ 

.702 

+     .36 

6 

- 

.33^ 

+ 

.^50 

+     .33 

7 

- 

.2*46 

H 

.627 

+  1.12 

8 

+ 

.024 

- 

.062 

+      .27 

9 

+ 

.137 

- 

.349 

+  2.53 

10 

+ 

.080 

- 

.205 

+  1.02 

11 

+ 

.007 

- 

.017 

+  1.15 

12 

- 

.232 

+ 

.590 

+     .60 

13 

-i 

.O65 

- 

.167 

+  1.43 

14 

- 

.28C 

+ 

.714 

+      .83 

15 

+ 

.107 

- 

.274 

+  2.16 

16 

- 

.368 

+ 

.937 

+      .65 
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TaojLe  VIII     lon.rrc,  Square;.:  Golution  ■   "ice's  l6  Points 

~|  10,000  Bf  1  Z  1   (a)     1(f) 

5a 

100  +2.17         -  lA.21         -  'i.od  -  15.66  1 

1  +0.54839      +1.66820  +7.21659 

10,000  5f  +  106.07      +  27.95  +  119.81  0 

-  93.05      -  23.71  -  102.55    +6.548 

+     13.02       +     4.2';-  +     17.26     +6.548       1 

-  1  -O.32565  -I.32565 

1  •  +  30.84  +     55.17 

-  6.04  -    26.12 

-  I.38  -      5.62 

[w]  +  23.42  +     23.43 


10,000  8f    Check     ^~      Check 


1  -0.32565   -I.32565   +1.66820  +7.21639 
10,000  6f                     -2.13248  -8.68087 

2  -0.32565   -1.32565   -0.46428  -1.46426 
Standard  Error  of  Single  Observations- 


» =  \/j™t  =yi#  =  ov88 

V  n  -  u    V  30 


_  6.548      _ 
Srf  "  13702      "     °^0i 

(+  2.17  Q   -  1^.21  (0.5^3)  -1=0 
8.148     Jt.- 
Qaa  =  -2717   =  3°755 

Standard  Errors  of  Unknowns  :  -     m  =  \x  J   Q 

m  

—^  =  0.8836  -J   3.755  =  1.712,   10,000  mf=0.8836\/o7o77=0. 24520 

m       =  171 o 2  m  m -  -  0.000024520 

a  1 
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Results  of  Least  Squares  Solutions  - 

10,000  6  f  =  -  0.32565      +  0.21*520 

5  f  =   -  0.000032565       +  0.0000^520 

f  =   0.00336T003 
o    

f  =  f     +  6  f  =  0.00333^38       +  0.000024520 

~  =  299.9     +     2.2 

5  a 

100   =  -  0.461+28  +  1.71  meters 

6  a   =      *  kS     +  171 

a0  =   6,378,388 
a  =  a  +  5  a   =   6,378,34-2  +  171  meters 

Here,  as  in  the  problem  of  section  7.1>  the  standard  errors 
of  the  corrections  are  almost  as  large  as  or  larger  than  the 
corrections  themselves.  But  the  fact  that  they,  as  well  as  the 
corrections  to  f  and  a  themselves,  are  relatively  small,  testifies  as 
to  the  value  of  the  method.  However,  it  must  be  realized  that  the 
standard  errors  reflect  only  the  precision  of  the  least  squares 
adjustment,  assuming  the  given  data  has  no  error,  which  of  course  is 
not  true. 
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7.3  Computation  of  the  Flattening  from  Gravity  Anomalies. 

Given  a  set  of  isostatic  gravity  anomalies  (Hayford, 
D  =  113. 7  tan.)  at  67  points  at  various  random  positions  scattered 
over  the  earth,  and  computed  on  the  basis  of  Helmert's  gravity 
formula  of  1901,  we  can  demonstrate  the  solution  of  a  problem  of 
correcting  the  gravity  formula  and  of  obtaining  a  value  for  the 
flattening  of  the  ellipsoid  (13),  p.  1&-19,    (1*0 .  The  method  is  that 
of  Chapter  5*  but  extended  to  JLnclude  the  longitude  term  of  the 
formula. 

Helmert's  Formulas - 

7  =  978.030  (1  +  0.005302  sin2cp  +  0.000007  sin22cp) 

for  which  the  p  term  yields  the  value  of  f  =  1/298.2.  See 
(13),  P.  78. 

The  general  form  of  the  gravity  formula,  including  the 

longitude  term  iss- 

2        2         2 
7  =  7„[1  +  P  sin  cp  +  e  sin  2cp  +  r  cos  cp  cos  2(\  -  X  )] 

a  o 

?  2  2 

=  7   [1  +  p  sin  cp+esin  2cp+r  cos  cp  (cos2\cos2X  +  sin  2\  sin  2\  )] 

Making  corrections,  we  have  7'  =  7  +  A7  , 

2        2  2 

where  A/  =  x1  +  y'sin  cp  +  z'cos  cp  cos  2\  +  u'cos  cp  sin  2\ 

x'  =  unknown  corr'n  to    7-,    in  mgals 

Hi 

y'  =    "      "    "      P     in  mgals  (•  1.022  x  lo"3) 
z'  =    "      "    "  r  cos  2\  in  mgals  (*  1.022  x  10~3) 


11  r  sin  2X  in  mgals  (♦  1.022  x  10"3) 
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Using  &f  -  £g  =   -  v,  we  have 

2  2  2 

x'    +  y'sin  cp  +  z'cos  cp  cos  2  X  +  u'cos  cp  sin  2  \  ■  /g=   -  v 

The  observation  equations  will  then  bes- 

ax'  +  b  y1  +cz'  +  d  u'  -  1  =  -  v 

where     a  =  1 

•  2 
b  =  sin  cp 

2 
c  =  cos  cp  cos  2  \ 

2 
d  =  cos  cp  sin  2  \ 

1  *  4g 

The  solution,  according  to  the  method  of  least  squares  will 

2 
be  Z  v  =  minimum. 
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Table  IX  -  Isostatic  Gravity  Anomalies  (Hayford,  D  =  113. 7  km.) 


Point 


9 


£% 


Point 


9 


4g 


1 

79?5 

11?5 

+ 

16 

34 

22?5 

36?5 

+  26 

2 

79-5 

17.5 

+ 

8 

35 

18.5 

38.5 

+  48 

3 

76.5 

15.5 

- 

26 

36 

14.5 

-    17.5 

+  44 

4 

71.5 

27.5 

+ 

20 

37 

13.5 

42.5 

+  ^3 

5 

65.5 

24.5 

- 

5 

38 

-  15.5 

38.5 

+  27 

6 

64.5 

11.5 

+ 

28 

39 

-  26.5 

15.5 

+  34 

7 

62.5 

17.5 

- 

20 

4o 

50.5 

120.5 

-     1 

8 

60.5 

22.5 

+ 

24 

41 

49.5 

94.5 

+    9 

9 

58.5 

38.5 

-H 

16 

42 

48.5 

116.5 

+  10 

10 

57.5 

4o.5 

+ 

28 

43 

47.5 

96.5 

+  19 

11 

55.5 

13-5 

+ 

10 

44 

47.5 

79.5 

+    6 

12 

54.5 

39.5 

- 

11 

45 

46.5 

123.5 

-    5 

13 

54.5 

55.5 

+ 

4 

46 

46.5 

92.5 

+  58 

14 

53.5 

29.5 

+ 

15 

47 

46.5 

60.5 

-    1 

15 

53.5 

50.5 

+ 

23 

48 

45.5 

102.5 

+  23 

16 

52.5 

25.5 

+ 

15 

49 

44.5 

110.5 

+  14 

17 

>l-5 

-      0.5 

+ 

9 

50 

44.5 

74.5 

+  26 

18 

51.5 

11.5 

+ 

36 

51 

43.5 

9^-5 

+  14 

19 

50.5 

10.5 

+ 

32 

52 

43.5 

74.5 

-  16 

20 

48„5 

16.5 

+ 

31 

53 

42.5 

71.5 

+    1 

21 

U7. 5 

39-5 

+ 

43 

54 

39.5 

120.5 

-  20 

22 

1+6.5 

7.5 

+ 

20 

55 

39.5 

84.5 

-  11 

23 

46.5 

34.5 

+ 

6 

56 

38.5 

78.5 

-    5 

24 

45.5 

10.5 

+ 

30 

57 

36.5 

87.5 

+  14 

25 

kk.  5 

-      1.5 

+ 

6 

58 

35.5 

105.5 

+  11 

26 

44.5 

39.5 

+ 

31 

59 

34.5 

79.5 

+  10 

27 

43.5 

42.5 

+ 

43 

60 

31.5 

94.5 

-    4 

28 

4o.5 

44.5 

+ 

74 

61 

56.5 

59.5 

+  55 

29 

39-5 

53.5 

+ 

15 

62 

51.5 

104.5 

-    2 

30 

36.5 

+      6.5 

+ 

14 

63 

46.5 

61.5 

+  14 

31 

^5.5 

14.5 

+ 

70 

64 

38.5 

70.5 

-  66 

32 

28.5 

33.5 

+ 

17 

65 

23.5 

79-5 

+  18 

33 

26.5 

33-5 

+ 

31 

66 
67 

22.5 
11.5 

77.5 
78.5 

+  27 

-  4o 
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Table  X  -  Error  Equations  -  Flattening  and  Correction 


a 

to  Gravity 

Formula  f rom 

Gravity 

Anomalies 

b 

c 

d 

1 

Point 

(x') 

(?) 

(z«) 

(u») 

=  -  £g 

1 

1 

.9669 

.0306 

.0130 

-  16 

2 

1 

.9669 

.0272 

.0190 

-  8 

3 

1 

.9^56 

.0467 

.0281 

+  26 

4 

1 

.8993 

■  0977 

.0825 

-  20 

5 

1 

.8281- 

.1128 

.1298 

+  5 

6 

1 

.8l47 

.1706 

.0724 

-  28 

7 

1 

.7868 

.1746 

.1223 

+  20 

8 

1 

.7576 

.1714 

.1714 

-  24 

9 

1 

.7269 

.06l4 

.2660 

-  16 

10 

1 

.7113 

.0452 

.2851 

-  28 

11 

1 

.6791 

.2858 

.1456 

-  10 

12 

1 

.6628 

.0643 

.3310 

+  11 

13 

1 

.6628 

-  .1209 

.3148 

-  4 

14 

1 

.6463 

.1822 

.3033 

-  15 

15 

1 

.6463 

-  .0062 

.3537 

-  23 

16 

1 

.6295 

.2332 

.2880 

-  15 

17 

1 

.6125 

.3874 

-  ,00#3 

-  9 

18 

1 

.6125 

.3567 

.1514 

-  36 

19 

1 

.595^ 

.3778 

.1450 

-  32 

20 

1 

.5610 

.3682 

.2391 

-  31 

21 

1 

.5^36 

.0871 

.4480 

-  43 

22 

1 

.5262 

.4577 

.1226 

-  20 

23 

1 

.5262 

.1698 

.4424 

-  6 

24 

1 

.5088 

.4586 

.1761 

-  30 

25 

1 

.4913 

.5081 

-  .0266 

-  6 

26 

1 

.4913 

.0971 

.4994 

-  31 

27 

1 

.^739 

.0459 

.5242 

-  43 

28 

1 

.4217 

.0101 

.5781 

-  74 

29 

1 

.4046 

-  .0726 

.5909 

-  15 

30 

1 

.3538 

.6297 

.1454 

-  14 

31 

1 

.3372 

.5796 

.3213 

-  70 

32 

1 

.2277 

.3017 

.7109 

-  17 

33 

1 

.1991 

.3129 

.7372 

-  31 
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Table  X     (cont'd) 


a 

b 

c 

d 

1 

Point 

(x«) 

(y) 

(*') 

(u') 

=  -  4 

3^ 

1 

.11*65 

.2U96 

.8163 

-  26 

35 

1 

.1007 

.2023 

.8763 

-  1*8 

36 

1 

.0627 

.7678 

-   .5376 

-  1+1* 

37 

1 

.051+5 

.0825 

.91+20 

-  1+3 

38 

1 

.071*** 

.2089 

.901+8 

-  27 

39 

1 

.1991 

.6865 

.k!2k 

-  31+ 

i*o 

1 

•  5951* 

-  .1962 

-  .3539 

+    1 

ki 

1 

.5782 

-  .1+165 

-   .0660 

-    9 

1+2 

1 

.5610 

-   .261*2 

-  .3506 

-  10 

1+3 

1 

.5^36 

-  .1*1*1*8 

-  .1027 

-  19 

kk 

1 

.5^36 

-  .1+261 

.1636 

-    6 

1+5 

1 

.5262 

-  .1852 

-   .1*362 

+    5 

k6 

1 

.5262 

-  .1+721 

-  .01*13 

-  58 

hi 

1 

.5262 

-    .21+1*1 

.1*062 

+    1 

1+8 

1 

.5088 

-    .1*1*52 

.2076 

-  23 

*+9 

1 

.1+913 

-  .381*0 

-  .3338 

-  11 

50 

1 

.1+913 

-  .1+361 

.2620 

-  26 

51 

1 

.1+739 

-  .5197 

-  .0823 

-  11* 

52 

1 

.1+739 

-  .1+511 

.2710 

+  16 

53 

1 

.1*561+ 

-  .1+31+1 

.3271 

-     1 

5h 

1 

.1+01*6 

-  .2886 

-   .5207 

+  20 

55 

1 

.1+01+6 

-  .581+1+ 

.1136 

+  11 

56 

1 

.3875 

-  .5638 

.2393 

+    5 

57 

1 

.3538 

-  .61*38 

.0561* 

- 11* 

58 

1 

•  3372 

-  .5681 

-   .31*13 

-  11 

59 

1 

.3208 

-  .63I+O 

.2U31* 

-  10 

bo 

1 

.2730 

-  .716C 

-  .1137 

+      4 

61 

1 

.6951* 

-  .li+Ti 

.2661+ 

-  55 

62 

1 

.6125 

-  .338c 

-  .1879 

+     2 

63 

1 

.5262 

-  .2581 

.3975 

-  11+ 

6k 

1 

.3875 

-  .1*759 

.3851* 

+  66 

65 

1 

.1590 

-  .7852 

.3011+ 

-  18 

66 

1 

.11*65 

-  .7736 

.3607 

-  27 

67 

1 

.0398 

-  .8839 

.3752 

+  1*0 
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Table  XI  -  Least  Squares  Solution  -  Flattening  and  Correction 

to  Gravity  Formula  from  Gravity  Anomalies 

Forward  Solution ;- 

x«     y»      z»      u'        1         L 

x'  +67     +  33-20     -     4.17     +  13.17     -      1033  -  923.80 

-  1     -0.4955     +O.0622     -O.1966     +  15.41Y9  +  13.7881 

y1  +  19.76     -     1.29     +     4.87     -  460.99  -  404.45 

-  16.45  +  2.07  -  6.53  +  511.85  +  457.74 

(1)  +  3.31  +  0.?8  -  1.66  +  50.86  +  53.29 

-  1    -0.2356  +0.5015  -  15.3656  -  16.0997 

z1  +  10.73  +  0.94  -  183.53  -  177.32 

-  0.26  +  0.82  -  64.25  -  57.46 

-  0.18  +  0.39  -  11.98  -  12.55 

(2)  +  10.29  +  2.15  -  259.76  -  247-33 

1    -O.2089  +  25.2439  +  24.0350 

u'  +  9.65  -  337.37  -  308.74 

-  2.59  +  203.09  +  181.62 

-  O.83  +  25.51  +  26.72 

-  0.45  +  54.26  +  51.67 

(3)  +  5.78  -  54.51  -  48.73 

-  1   +  9.4303  +  8.4308 

1  +52,563  +50,548.11 

-15,926.69  -14,243.06 

-  781.49  -   818.83 

-  6557.36  -  6243.32 

514.07  -       459.56 

[w]  +28,783.39  +28,783.34 


Rack  Solution: 


Check     y*      Check     x*      Check 


1  +9.4308  +8.4308  +25.2439  +24.0350  -15.3656  -16.0997  +15.4179  +13.7881 

u*  -  1.9701  -  1.7612  +  4.7295  +  4.2280  -  1.8541  -  I.6575 

z»  -  5.4833  -  5.2477  ♦  1.4476  +  1.3854 

y'  +  f.9872  +  8.4827 

Z  +9.4308  +8.4308  +23.2738  +22.2738  -16.1194  -17.1194  +22.9986  +21.9987 


75 


N 


IP 
I 

U 


o 
■p 


a 
-p 


M 


H 

O 

o 

ON 

CO 

O 

O  OJ 

• 

o 

i 

ON 

$ 

OJ 

• 

o 

1 

H 

VO 

IfN 

CO 

O  O  OJ 

0 

o 

1 

VO 
IfN 

on 

OJ 

9 

o 

1 

o 

o 

IfN  OJ 
r-1  ON 
O-* 
IfN  O 

d  d 
+  + 

o 

IfN 
IfN 

• 

o 
+ 

iH 

IfN 
0-=fr 

IfN 

&> 
ON 

OJ  t— 
OJVO 
VO   rH 

OOH 

ON 

CO 
t- 
H 

o 

vo  ifN-4- 

VO  en  O- 

ON-d-  on 

H  OJ  O 

IfN 
OJ 

3 

• 

o 

1 

• 

O 
I 

o  o 

+  + 

• 

O 

+ 

o 

O 

« 

0          o 

o  o 

1     1 

• 

o 

g 

H  ON 

oo    • 
H  O 
+     l 

VO  H 

vo  o 

•   IfN 

o 
I  + 

ON 

•  OJ 
OJ      • 

o 
+    1 

• 

IfN 

+ 

C-  OJ 

J*    . 

o 

1  + 

co  iR 

C—  CO 

A  ™ 

o    . 
o 

+   1 

ON 

OJ 

O 
r4 

+     I 

O  ITN 

OJ  IfN 

.  ON 

CO    • 

O 

+     i 

H 
PO 

.  H 

CO 

+    1 

C-  H 

X 


N 


OJ 


PO 


76 


Solution  of  Correlate  Numbers :- 

5.78  Q   =  -  0.14-825 
xu 

ft   =  -  O.08W 

10.29    0x2    =    +    2'1^>     ("     .08478)    =    0.1789 

Q^z  =  (0.1789  +  0.1823)    *  10.29  =   0.03510 

3.31  Q  +  O.78  (.03510)  -  1.66  (-  .08478)  =  -  0.4955 

Q       =  (-  0.4955  -  0.0274  -  0.1407)   *  3.31  -   -  0.20048 

67  Q       +  33.20  (-  .20048)   -  4.17   (.03510)  +  13.17  (-  .08478)  =  1 

0      -  (1  +  6.6559  +  0.1464  +  1.1166)  ♦  67  =  0.13312 

5,78  %i  =  °'5507 
Q^  =  0.09528 

10.29  Q      +  2.15  (.09528)  =  -  0.2356 

0      =  (-  0.2356  -  0.2049)  *  10.29  =  -  0.04281 

3.31  Q      +  0.78  (-  .04281)  -  1.66  (.09528)  =  1 

ft      =  (1  +  0.0334  +  0.1582)  1  3.31  =  O.36000 

5.78  Q       =   -  0.2089 
zu 

Q       =    -  0.03614 
zu 

10.29  Q       +  2.15   (-  .03614)  =  1 

qzz  =  (1  +  0.0777)  *  10.29  =  0.10473 
5.78  ft^  =  1 

0       =  0.17301 

A1U 
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Standard  Error  of  an  observation  of  unit  weight:- 

u   =^n^]u  =  ^  W^T    =  J  ^6'87  =  ±  21.37  fflgal 
Standard  Errors  of  the  unknowns  after  adjustment :- 


m' 


l'x  =  ^>7"~^x~'  etc- 
m  *  =  21.37  V  0.1331  =  +  7.78  (mgal. ) 

.A.  mm 

m  ■  =  21.37  n/  0.3600  =  +  12.82  (mgal.  *  1.022  x  lo"3) 
1  =  21.37  v/  0.1047  =  +  6.92  (mgal.  *  1.022  x  10"3) 
•   -  21.37    n/  0.1730     =    %+     8.89  (mgal.    f  1.022  x  lo"3) 


m 
z 


m 
u 


The  values  for  the  unknowns,  solved  in  Table  XI,  and  their 
standard  errors,  are  here  converted  to  gals.:- 

x  =  x1  x  lo"3  =  +  0.0230   +  0.0078  gals, 

y  =  y1  x  (1.022)  x  10"  =  -  0.000016  +  0.000013  gals, 

z  =  z»  x  (1.022)  x  10"  =  +  0.0000238+  0.0000071  gals, 

u  =  u1  x  (1.022)  x  10"  =  +  O.OOOOO96+  0.0000091  gals. 

The  components  of  the  longitude  term  are:- 

r2   =  z2  +  u2  =  O.OOO6586  x  10"  ,  r  =  0.0000257 


tan2^o=  ?■  =1OT  =  °-Uo52 


2  \     =  22  ,   X     =   11 
o     '    o 


Applying  the  computed  corrections : - 

7E      =  978.03C  +  x  =  978.053   gals, 
p    -  0.005302  +  y  =  0.005286 
e    =  0.000007      (unchanged) 
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and  adding  the  longitude  term:- 

0.000026  cos2cp  cos  2  (X  -  11°) 

we  obtain  the  corrected  gravity  formula: - 

7   =  978.053  [1  +  0.005286  sin2cp  +  0.000007  sin2  2  cp 
+  0.000026  cos2cp  cos  2  (X   -  11°)]. 

Next,  the  flattening  f  is  solved  from  the  corrected  value  of  p. 
At  this  point  it  should  be  pointed  out  that  if  the  principal  problem 
is  to  find  the  flattening,  it'would  not  be  necessary  to  include 
solution  of  the  longitude  term.  With  slight  difference  the  same 
result  for  the  correction  to  (3  is  obtained  when  only  two  unknowns  exist 
in  the  normal  equations.  Witness,  using  the  first  reduced  normal 
equation  in  Table  XI:- 

3.31  y'  +  50.86  =  0 

y1   =   -  50.86   *  3.31  =   -  l%k 
x       =   -  15.^  x  1.022  x  10"     =   -  0.000016 
0       =  0.005302   -  0.000016  =  0.005286 

Now  using  Clairaut's  Formula,    (5-1)*  we  obtain  f  from  the  new 
value  of  0. 

0o  0052855     =     |m-f-i£mf 

oTa  2  7T  -1 

m      -  77  *      =  S0T0T  sec 

E 

/  27T  >,2     637,^20,000  cm. 
m       ~    ^86T6V  970.05       galT 

-  3?.  W  x  637.82  25,179.858  nn^7H 

"  7^.2  x  978.05  7,26l,12U     "     °-°°3^7« 
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0.0052855  =  |  (0.003^678)  -  f  [1  +  i|  (0.003^78)] 
=  0.008670    -    f  (1  +  o.ooi+2109) 

O.OO8670   -  0.005286 
1.0042109 

=  g;gSai09    =     °-°°33698    +    O.OOOOO83 

l/f        =  296.75    +     0.7^ 

(3  =  0.0052855  +  0.0000130 
:.  mf        =+ m   (0.0000131)  ■  +  O.OOOOO83) 


Of  course,  the  result  reflects  the  fact  that  only  a  few  values 
of  £g  were  used,  and  at  scattered  locations.  The  solution  by 
Heiskanen,  referred  to  in  Chapter  5>  shows  that  the  use  of  a  great  many 
points  results  in  a  value  of  f  equal  to  that  obtained  from  Helmert' s 
original  formula.  The  results  in  the  computation  given  here  are  sur- 
prising considering  the  fact  that  the  points  chosen  were  completely 
random  with  no  thought  given  to  even  approximately  even  distribution. 

7.4  Computation  of  Gravimetric  Deflections 

7.^1  Procedure  and  General  Information. 

The  procedure  used  for  computation  of  the  gravimetric  deflect- 
ions at  the  two  points  chosen,  is  that  involving  the  use  of  the  Rice 
Circle -Rings.  The  extent  of  readily  available  gravity  material 
dictated  that  the  area  covered  in  both  cases  was  to  be  wtthic  the 
circle-ring  no.  h$,   having  a  radius  of  399*0  km. 
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A  template  was  prepared  to  the  scale  of  the  available  anomaly 

o 

maps,  having  circles  of  radii  listed  in  Table  I  and  with  10  sectors. 

The  usual  manner  of  orienting  the  template,  which  is  to  have  zero 
azimuth  bisect  one  of  the  sectors  making  mean  azimuths  of  all  sectors 
multiples  of  10  ,  was  not  followed  here  because  the  writer' s  method  of 
assuring  proper  sign  determination  dictated  that  each  compartment  have 
its  own  quadrant  designation.  Therefore  the  azimuth  lines  were 
numbered  in  multiples  of  10  ,  and  the  mean  sector  azimuths  are  multi- 
ples of  5  • 

For  the  point  in  the  northern  hemisphere,  zero  azimuth  is  nortij, 
and  for  the  one  in  the  southern  hemisphere  it  is  south.  Following  the 
usual  practice  for  points  in  Europe,  \   and  t\   were  considered  positive 
eastward,  so  for  both  points,  the  azimuths  were  measured  eastward  from 
zeroj  thus  they  are  measured  clockwise  in  the  northern  hemisphere,  and 
counterclockwise  in  the  southern  hemisphere.  This  will  preserve  signs 
automatically,  using  the  same  procedures  for  both.  Figure  12  illu- 
strates the  azimuth  convention  used. 

Remembering  that  the  plumb  line  (below  the  surface)  is  pulled 
toward  the  mass  surplus  and  thus  in  the  direction  of  the  more  positive 
of  the  anomalies  on  opposite  sides  of  the  computation  point,  we  can 
devise  a  system  of  assuring  correct  sign.  Figure  13  illustrates  this. 
|  is  considered  positive  toward  the  nearest  pole,  t\   is  considered 
positive  eastward.  Table  XIII  gives  criteria  to  be  followed. 
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Fig.    12 


Azimuth  Convention  for  Computing 
Gravimetric  Deflections  for  the 
Points  Computed  in  this  Thesis. 
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Fig.   13 
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Table  XIII  -  Criteria  for  Determining  Effective  Azimuth 

and  Sign  of  Deflection 


Quadrants 

Anoma.ly  Gradient 

Effective 
Azimuth 

cosine 
and 

g 

sine 
and 

n 

I  and  III 

(Agj  -  4gXII)  <  o 

*i 

>   0 

>   0 

(^  -  Agi;[I)  >  0 

*in 

<  0 

<  0 

II  and  IV 

(^n-AgJV)  <o 

*ii 

<  0 

>  0 

(AgZI  -  AgIV)  >  0 

^v 

>  0 

<  0 

. 

i 

Computation  proceeds,  "beginning  in  the  first  quadrant  with  the 
effect  of  that  quadrant  and  the  one  opposite,  the  third.  After  com- 
pleting the  first  and  third,  the  second  quadrant  is  begun,  determining 
the  effect  of  that  quadrant  and  the  one  opposite  it,  the  fourth.  Thus 
the  entire  3^0  in  azimuth  are  covered. 

In  order  to  obtain  data  for  a  graph  of  the  effects  of  areas  of 
different  radii  on  the  deflection  of  a  point,  the  computation  was  made 
grouping  the  effects  of  four  zones  together  from  the  point  outward. 
Then  accumulative  effects  were  computed  and  plotted  for  both  £  and  rj. 

The  inner  circle,  for  one  case  5*125  km.  in  radius,  and  for  the 
other  10.15  km.,  was  computed  by  the  one-gradient  method.   Insufficient 
local  gravity  data  at  the  points  dictated  that  the  three -gradient 
method  would  yield  results  of  meaningless  accuracy. 

The  two  points  computed  were  Laiska,  Finland  (6k     03'  N, 
28°  50'  E),  and  Wonderfontein,  Union  of  South  Africa  (25°  WJ1  S, 


8>* 


o      \ 
29  53'  E).  The  points  were  chosen  to  contribute  data  for  the  com- 
putation of  the  equatorial  radius  from  the  arc  of  the  30th  meridian. 
In  addition,  their  locations  had  to  be  such  that  the  anomaly  maps 
available  to  the  writer  would  provide  data  to  a  large  enough  radius 
around  each  point  to  make  computation  by  the  circle -ring  method 
feasible.  The  result  was  that  both  the  points  turned  out  to  be  some- 
what removed  from  areas  of  concentrated  gravity  stations.  However, 
point -anomaly  values  obtained  from  files  of  the  Institute  of  Geodesy, 
Fhotogrammetry,  and  Cartography,  The  Ohio  State  University,  for  the 
general  area  around  the  points,  allowed  the  plotting  of  iso-anomaly 
contour  maps  to  the  scale  1"  =  6  km.  (l  :  2*4-0,000)  for  both  points,  and 
the  writer  proceeded  with  the  computation  undaunted  by  the  lack  of 
enough  gravity  data. 

The  small  scale  maps  used  for  the  bulk  of  the  computation  were 
the  following: - 

Laiska  :  Free -air  anomaly  map  of  Finland  on  the  scale 

1  :  1,000,000,  covering  the  area  21  -  35  East  longitude 
(extending  into  Russia)  and  60  -  69  North  latitude. 
Contours  for  the  Gulf  of  Bothnia  were  sketched  in  from  a 
smaller  scale  map  of  that  area.  The  maps  used  were 
prepared  from  gravity  determinations  of  the  Finnish 
Geodetic  Institute  and  were  computed  on  the  International 
Gravity  Formula.  Contour  interval  was  5  mgal. 

Wonderfontein  :  Bouguer  anomaly  lithographed  map  sheets  of 
Union  of  South  Africa,  Northeast  and  Southeast  sheets,  on 
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scale  1  t  1,000,000,  10  mgal  contours;  total  coverage  of 
the  two  sheets  s  22°  00'  to  3^°  30'  S,  23°  30'  to  33°  50' 
E.  These  maps  were  published  in  South  Africa. 

Topographic  maps  published  as  World  Aeronautical 
Charts  by  U.S.  Air  Force  (ACIC),  scale  1  :  1,000,000; 
contour  interval  is  1000  ft.,  but  adequate  point  elevations 
to  the  nearest  foot  are  given  for  the  use  to  which  the  maps 
were  put.  Total  coverage  of  sheet  numbers  1275*  127&,  1299* 
1300,  1397*  and  1398  is  about  the  same  as  that  of  the 
Bouguer  aaomaly  maps. 

Since  the  South  African  point  was  estimated  from  a  Bouguer 
anomaly  map,  it  was  necessary  to  use  the  topographic  maps  also,  for  the 
purpose  of  obtaining  mean  elevations  for  each  compartment  of  the  same 
Rice  ring  template  used  for  the  anomalies.  Thus  armed  with  the 
elevations,  the  Bouguer  correction  could  be  computed  and  subtracted 
from  the  Bouguer  anomaly  to  obtain  the  Free  Air  anomaly  for  the  com- 
partment. 

<figf  =  AgB  -  g'   =  Free  Air  anomaly 

where  z^^   =  Bouguer  anomaly 
a 

g1   =  Bouguer  correction 
B 

3  P 


ir<2§>h 


m 


p   =  density  of  crust  assumed  for  area 
=  2.67 


p   =  mean  density  of  earth 
m 


=  5.576 
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/.  g'B  =  -|  |f|^  (0.09^06)  h  mgals     (h  in  feet) 
=  -  0.03378  h  mgals. 


7.^2  Computation  for  Laiska,  Finland. 

The  computations  for  the  circle-ring  zones,  proceeding  from  the 
outside  zone,  no.  k8,   and  continuing  inward  toward  the  computation  point 
to  zone  no.  27,   are  included  in  Table  XIV. 

The  inner  circle,  r  =  10.15  km«>  is  computed  below.  The 
effects  of  both  the  inner  circle  and  the  circle-ring  zones  are 
accumulated  in  Table  XV.  Curves  of  the  effect  on  deflections  of  in — 
creasing  gravity  field  radius  are  shown  in  figure  15 . 

tN 
+2.0 


-0.3 


+--I U4-.S 


+ZZ 


Fig.  Ik     -  Inner  Circle  -  Laiska. 


Using  the  values  shown  in  fig.  Ik,   we  obtain: 


d|"  =  0.0525  (^gs  -  4gN) 

=  0.0525  (2.2  -  2.0)  =  0.0525  (0.2)  =  +  07 010 

dT)"  =  0.0525  (£gw  -  £gE) 

=  0.0525  (-  0.3  -  ^.8)  =  0.0525  (-  5.1)  =  -  05268 
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Table  XIV  -  Computation  of  Gravimetric  Deflections 

at  Laiska,  Finland 

£g  (mgal) 


Quadrant 

Sector 

Zone 

I 
0-10 

III 
180-19C 

• 

10   ■ 

c 

-  20 

III 
190-200 

I 
20  -  30 

III 
200-210 

w 

-     10 

-    30 

15 

-     30 

-     20 

-     20 

hi 

-     20 

-      5 

- 

25 

-     15 

-     10 

0 

L  £g 

-    30 

-    35 

- 

1+0 

-     ^5 

-     30 

-     20 

I  -   III 

+      5 

+      5 

-     10 

d0  @  Ol'OOl 

o'.'oo5 

0"005 

0V010 

Eff.   A 

185 

195 

25 

cos  A 

-  .9962 

-  .9659 

+ 

.9063 

Pin  A 

-  .0872 

-  .2588 

+ 

.1+226 

d| 

-   0'.'005 

-  0"005 

+ 

0"  009 

dt) 

0V000 

-  O'.'OOl 

+ 

o'.'ooi+ 

1+6 

-    15 

+  15 

- 

10 

+     20 

0 

+     20 

^5 

0 

+  25 

+ 

5 

+     20 

+      5 

+     20 

kh 

+     10 

+  30 

+ 

10 

+     15 

+     10 

+     10 

43 

+     10 

+     20 

+ 

10 

+     15 

+  15 

+     10 

z 

+      5 

+     90 

+ 

15 

+     70 

+  30 

+     60 

I  -  III 

-    85 

- 

55 

-    30 

d0 

O.085 

0. 

055 

0.030 

A 

5 

15 

25 

cos 

+  .9962 

+  .9659 

+ 

.9063 

sin 

+  .0872 

+  .2588 

+ 

.1+226 

d£ 

+     .085 

+     . 

053 

+ 

.027 

dri 

+     .007 

+     . 

Oll+ 

+ 

.013 

1+2 

+    15 

+     15 

+ 

15 

+     15 

+     18 

+       7 

1+1 

+    18 

+     11 

+ 

20 

+    13 

+     22 

+     10 

1+0 

+    18 

+     10 

+ 

20 

+     10 

+     25 

+    11+ 

39 

+    15 

+       8 

+ 

20 

+      5 

+     2l+ 

+      7 

z 

+    66 

+     kk 

+ 

75 

+    43 

+     89 

+     38 

I  -  III 

+     22 

+    32 

+    51 

de 

0.022 

0.032 

0.051 

A 

185 

195 

205 

cos 

-   .9962 

-  .9659 

-  .9063 

sin 

-   .0872 

-  .2588 

-     1+226 

d| 

-      .022 

-     .031 

-     .01+6 

dn 

-      .002 

-     .008 

-     .022 
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Table  XIV     (cont'd) 


I 

III 

I 

III 

I 

III 

30  - 

-  1+0 

210-220 

1+0  - 

50 

220-230 

50  - 

60 

230-21+0 

1+8 

+ 

5 

+     10 

+ 

20 

0 

+ 

35 

-     10 

hi 

+ 

10 

+     10 

+ 

30 

-      5 

+ 

^5 

-     10 

z 

+ 

15 

+     20 

+ 

50 

-      5 

+ 

80 

-     20 

I  -  III 

- 

5 

+    55 

+   100 

de 

005 

.055 

.100 

A 

35 

225 

235 

cos 

+  .8192 

* 

-  .7071 

-  .5736 

sin 

+  -5736 

-  .7071 

-  .8192 

ag 

+    . 

.00*+ 

-     .039 

-   .057 

dTj 

+    , 

003 

-     .039 

-    .082 

k6 

+ 

5 

+      5 

+ 

20 

+      5 

+ 

25 

+     10 

^ 

+ 

5 

+      5 

+ 

10 

+    10 

+ 

15 

+     5 

kk 

+ 

10 

+     10 

+ 

10 

0 

+ 

10 

0 

h3 

+ 

15 

+      5 

+ 

10 

-      5 

+ 

10 

0 

Z 

+ 

35 

+    25 

+ 

50 

+     10 

+ 

60 

+    15 

I  -  III 

+    10 

+     1+0 

+    1+5 

d0 

.010 

.Ol+O 

.01+5 

A 

215 

225 

235 

cos 

-  .8192 

-  .7071 

-  .5736 

sin 

-  .5736 

-  .7071 

-  .819? 

d| 

-    .008 

-    .028 

-    .026 

cLrj 

-    .006 

-    .028 

-     .037 

1+2 

+ 

18 

+      6 

+■ 

17 

+      5 

+ 

11+ 

+      5 

in 

+ 

22 

+     12 

+ 

22 

+     12 

+ 

17 

+     10 

1+0 

+ 

27 

+     15 

+ 

26 

+     15 

+ 

20 

+      5 

39 

+ 

3^ 

+     13 

+ 

28 

+     13 

+ 

20 

+      6 

z 

+ 

101 

+    1+6 

+ 

93 

+    1+5 

+ 

71 

+    26 

I  -  III 

+    55 

+    1+8 

+    1+5 

d0 

.055 

.01+8 

.01+5 

A 

215 

225 

235 

cos 

-  .8192 

-  .7071 

-  .5736 

sin 

-  .5736 

-  .7071 

-  .8192 

d| 

-     .0^5 

-     .03^ 

-     .026 

dt) 

-     .032 

-     .03^ 

-     .037 
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Table  XIV     (cont'd) 


I 

III 

I 

III 

I 

III 

6o  - 

70 

240-250 

70  - 

80 

250-260 

80  - 

90 

260-270 

48 

+ 

^5 

-  15 

+ 

45 

-  4o 

+ 

35 

-  4o 

47 

+ 

4o 

-  15 

+ 

35 

-  25 

+ 

30 

-  35 

Z 

+ 

85 

-  30 

+ 

80 

-  65 

+ 

65 

-  75 

I  -  III 

+  115 

+  145 

+  l4o 

&G 

.115 

.145 

.l4o 

A 

245 

255 

265 

cos 

-  .4226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

*t 

-  .049 

-  .038 

-   .Oiii 

dTl 

-  .104 

-  .140 

-   .139 

46 

+ 

15 

-  10 

+ 

15 

-  20 

+ 

20 

-  25 

45 

+ 

10 

0 

+ 

10 

-  15 

+ 

15 

-  25 

44 

+ 

10 

+   5 

+ 

10 

-  10 

+ 

20 

-  10 

43 

+ 

10 

-   5 

+ 

15 

-   5 

+ 

20 

-   5 

z 

+ 

45 

-  10 

+ 

50 

-  50 

+ 

75 

-  65 

I  -  III 

+  55 

+  100 

+  l4o 

de 

.055 

.100 

.140 

A 

245 

255 

265 

cos 

-  .4226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

d£ 

-  .023 

-  .026 

-   .012 

dri 

-  .050 

-   .097 

-   .139 

42 

+ 

13 

+   5 

+ 

17 

0 

+ 

15 

-   3 

ki 

+ 

12 

+   5 

+ 

10 

+   3 

+ 

5 

0 

40 

+ 

10 

0 

+ 

5 

0 

+ 

4 

0 

39 

+ 

10 

0 

+ 

6 

0 

+ 

7 

0 

z 

+ 

45 

+  10 

+ 

38 

+   3 

+ 

31 

-   3 

i  -  in 

+  35 

+  35 

+  34 

de 

.035 

.035 

.034 

A 

245 

255 

265 

cos 

-  .4226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

as 

-  .015 

-  .009 

-  .003 

dT) 

-  .032 

-  .034 

-  .034 

90 


Table  XIV     (cont'd) 


ii 

IV 

II 

IV 

II 

IV 

90  - 

100 

270-280 

100- 

110 

280-290 

110- 

120 

290-300 

k8 

+ 

30 

-  15 

+ 

15 

-     20 

+ 

10 

-     20 

47 

+ 

20 

-  30 

+ 

10 

-    30 

+ 

5 

-    35 

t 

+ 

50 

-  45 

+ 

25 

-     50 

+ 

15 

-     j>1> 

II  -  IV 

+    95 

+    75 

+     70 

de 

.095 

.075 

.070 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.4226 

sin 

- 

.9962 

- 

.9659 

- 

.9063 

d| 

+ 

.008 

+ 

.019 

+ 

.030 

dri 

- 

.095 

- 

.072 

- 

.063 

k6 

+ 

15 

-    30 

+ 

5 

-    35 

0 

-     25 

45 

+ 

10 

-    25 

0 

-    25 

- 

5 

-     25 

44 

+ 

15 

-      5 

0 

-      5 

- 

10 

-     20 

43 

+ 

10 

+      5 

- 

5 

+      5 

- 

5 

-     15 

Z 

j 

50 

-    55 

0 

-    60 

• 

20 

-     85 

II  -  IV 

+  105 

+    60 

+     65 

de 

.105 

.060 

.065 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.4226 

sin 

- 

.9962 

- 

.9659 

- 

.9063 

&S 

+ 

.009 

+ 

.016 

+ 

.027 

dTj 

- 

.105 

- 

.058 

- 

.059 

42 

0 

2 

0 

-      7 

0 

-      5 

41 

+ 

7 

2 

+ 

12 

-      7 

+ 

12 

-      8 

40 

+ 

15 

2 

+ 

20 

-      6 

+ 

18 

-     10 

*9 

+ 

15 

1 

+ 

20 

-      7 

+ 

20 

-      8 

E 

+ 

37 

-       7 

+ 

52 

-    27 

+ 

50 

-    31 

II  -  IV 

+     44 

+    79 

+    81 

de 

.044 

.079 

.081 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.4226 

sin 

- 

.9962 

- 

.9659 

- 

.9063 

at 

+ 

.004 

+ 

.020 

+ 

.034 

dti 

- 

.044 

- 

.076 

- 

.073 
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Table  XIV     (cont'd) 


II 

IV 

II 

IV 

II 

IV 

120-130 

300-310 

130- 

140 

310-320 

140-150 

320-330 

48 

+ 

15 

-  15 

+ 

15 

-  20 

+ 

10 

-  10 

47 

+ 

10 

-  20 

+ 

15 

-  15 

+ 

15 

-  10 

E 

+ 

25 

-  35 

+ 

30 

-  35 

+ 

25 

-  20 

II  -  IV 

+  6o 

+  65 

+  45 

6.6 

.060 

.065 

.045 

A 

305 

315 

325 

cos 

+ 

•  5736" 

+ 

.7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

dg 

+ 

.034 

+ 

.046 

+ 

.037 

dri 

- 

.049 

- 

.046 

- 

.026 

46 

+ 

5 

-  25 

+ 

10 

-  15 

+ 

15 

-  15 

45 

- 

5 

-  25 

+ 

5 

-  15 

+ 

15 

-  10 

44 

- 

5 

-  20 

0 

-   5 

+ 

10 

-  10 

43 

0 

-  10 

0 

-   5 

+ 

10 

-  10 

£ 

_ 

5 

-  80 

+ 

15 

-  40 

+ 

45 

-  45 

II  -  IV 

+  75 

+  55 

+  90 

de 

.075 

.055 

.090 

A 

305 

315 

325 

cos 

+ 

•  5736 

+ 

.7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

d| 

+ 

.043 

+ 

.039 

+ 

.074 

dn 

- 

.061 

- 

.039 

- 

.052 

42 

+ 

7 

-   5 

+ 

5 

-   5 

+ 

8 

-   15 

41 

+ 

13 

-   8 

+ 

13 

-  10 

+ 

10 

-  20 

40 

+ 

17 

-  14 

+ 

16 

-  14 

+ 

10 

-  10 

39 

+ 

20 

-   8 

+ 

18 

-   5 

+ 

10 

-   8 

z 

+ 

57 

-  35 

+ 

52 

-  34 

+ 

38 

-  53 

II  -  IV 

+  92 

+  86 

+  91 

d0 

.092 

.086 

.091 

A 

305 

315 

325 

cos 

+ 

.5736 

+ 

.7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

d| 

+ 

.053 

+ 

.061 

+ 

.075 

dT] 

- 

.075 

- 

.061 

- 

.052 
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Table  XIV     (cont'd) 


II 

IV 

II 

IV 

II 

IV 

150- 

■16c 

330-3^0 

160- 

170 

340-350 

170- 

180 

350-360 

48 

0 

0 

. 

5 

-  10 

. 

10 

-  20 

47 

+ 

10 

-   5 

+ 

10 

-  15 

+ 

10 

-  15 

L 

+ 

10 

-   5 

+ 

5 

-  25 

0 

-  35 

II  -  IV 

+  15 

+  30 

+  35 

de 

.015 

.010 

.035 

A 

335 

345 

355 

cos 

+ 

.9063 

+ 

.9659 

+ 

.9962 

sin 

- 

.4226 

- 

.2588 

- 

.0872 

as 

+ 

.014 

+ 

.029 

+ 

.035 

at] 

- 

.006 

- 

.008 

- 

.003 

46 

+ 

15 

-   5 

+ 

15 

-  10 

+ 

25 

-   5 

45 

+ 

15 

-  10 

+ 

15 

-  10 

+ 

25 

0 

kk 

+ 

10 

-  20 

+ 

10 

-  10 

+ 

15 

+  10 

43 

+ 

5 

-  20 

+ 

5 

0 

+ 

10 

+  15 

L 

+ 

^5 

-  55 

+ 

45 

-  30 

+ 

75 

+  20 

II  -  IV 

+  100 

+  7^ 

+  55 

d0 

.100 

.075 

.055 

A 

335 

345 

355 

cos 

+ 

.9063 

+ 

.9659 

+ 

.9962 

sin 

- 

.4226 

- 

.2588 

- 

.0872 

d£ 

+ 

.091 

+ 

.072 

+ 

.055 

dt) 

- 

.042 

- 

.019 

- 

.005 

42 

+ 

8 

-  20 

+ 

7 

-  10 

+ 

9 

+  12 

4l 

+ 

10 

-  15 

+ 

10 

-   13 

+ 

12 

+   5 

40 

+ 

5 

-   8 

+ 

10 

-   5 

+ 

14 

+   4 

39 

+ 

4 

-   8 

+ 

3 

2 

+ 

10 

+   5 

z 

+ 

27 

-  51 

+ 

30 

-  30 

+ 

45 

+  26 

II  -  IV 

+  78 

+  60 

+  19 

de 

.078 

.060 

.019 

A 

335 

345 

355 

cos 

+ 

.9063 

+ 

.9659 

+ 

.9962 

sin 

- 

.4226 

- 

.2588 

- 

.0872 

as 

+ 

.071 

+ 

.058 

+ 

.019 

dT) 

- 

.033 

- 

.016 

- 

.002 

93 


+  19 

+  97 

+  18 

+  81 

+  23 

+  51 

+  79 

+  58 

.053 

.079 

.058 

185 

195 

205 

-  .9962 

-  .9659 

-  .9063 

-  .0872 

-  .2588 

-  .4226 

-  .053 

-  .076 

-   .053 

-  .005 

-  .020 

-   .024 

Table  XIV     (cont'd) 


I  III  I  III  I  III 

0-10       180-190       10  -  20       190-200       20  -  30       200-210 


38  +17  +9           +29           +6           +     2k           +3 

37  +16  +9          +25          +7          +20          +5 

36  +18  +1          +    2k         +5          +19          +9 

35  +21  0+19+1+18+6 

Z  +72 

I  -  III 

ae 

A 

cos 
sin 
ag 

at] 

3k  +16  0+19  3+18+2 

33  +12  0+16  0+17+2 

32  +10  +2          +13  +1          +16          +3 

31  +8  +      k          +12  +3          +16          +2 

E 

I  -  III 

de 

A 

cos 
sin 
dg 

dn 

30 

29 

28 

27 

z 

I  -  III 

de 

A 

cos 

sin 

ag 

dr) 


+ 

k6 

+   6 
+  ko 

.0+0 

185 

-  .9962 

-  .0872 

-   .0+0 

-  .003 

+ 

60 

+   + 

+  56 

.056 

195 

-  .9659 

-  .2588 

-  .054 

-  .01+ 

+ 

67 

+   9 

+  58 

.058 

205 

-  .9063 

-  .4226 

-  .053 

-  .025 

+ 

8 

+   5 

+ 

10 

+  k 

+ 

15 

+   3 

+ 

7 

+   5 

+ 

9 

+  + 

+ 

11 

+   3 

+ 

6 

+     k 

+ 

7 

+   3 

+ 

9 

+   3 

+ 

5 

+     k 

+ 

6 

+   3 

+ 

7 

+   3 

+ 

26 

+  18 

+   8 

.008 

185 

-  .9962 

-  .0872 

-  .008 

-  .001 

+ 

32 

+  1+ 

+  18 

.018 

195 

-  .9659 

-  .2588 

-  .017 

-  .005 

+ 

k2 

+  12 

+  30 

.030 

-  .9063 

-  .4226 

-  .027 

-  .013 

94 


Table  XIV     (cont'd) 


i 

III 

I 

III 

I 

III 

30  - 

40 

210-220 

40  - 

50 

220-230 

50  - 

60 

230-240 

38 

+ 

29 

+   7 

+ 

25 

+   8 

+ 

16 

+   5 

37 

+ 

22 

+   5 

+ 

18 

+   3 

+ 

15 

+   1 

36 

+ 

18 

+   5 

+ 

17 

+   2 

+ 

15 

+   1 

35 

+ 

17 

+   6 

+ 

15 

+   5 

+ 

14 

+   5 

E 

+ 

86 

+  23 

+ 

75 

+  18 

+ 

60 

+  12 

I  -  III 

+  63 

+  57 

+  48 

de 

.063 

.057 

.048 

A 

215 

225 

235 

cos 

-  .8192 

-  .7071 

-  .5736 

sin 

-  .5736 

-  .7071 

-  .8192 

&t 

-  .052 

-  .040 

-  .028 

dri 

-  .036 

-  .040 

-  .039 

34 

+ 

15 

+   6 

+ 

14 

+   7 

+ 

13 

+   8 

33 

+ 

15 

+   6 

+ 

14 

+   8 

+ 

14 

+  10 

32 

+ 

17 

+   4 

+ 

18 

+   l 

+ 

19 

+   4 

31 

+ 

19 

+   2 

+ 

20 

+   1 

+ 

20 

0 

S 

+ 

66 

+  18 

+ 

66 

+  17 

+ 

66 

+  22 

I  -  III 

+  48 

+  49 

+  44 

de 

.048 

.049 

.044 

A 

215 

225 

235 

cos 

-  .8192 

-  .7071 

-  .5736 

sin 

-  .5736 

-  .7071 

-  .8192 

at 

-  .039 

-  .035 

-  .025 

dT) 

-  .028 

-  .035 

-  .036 

30 

+ 

16 

+   2 

+ 

16 

+   1 

+ 

17 

+   1 

29 

+ 

12 

+   2 

+ 

12 

+   2 

+ 

13 

+   2 

28 

+ 

9 

+   3 

+ 

9 

+   2 

+ 

8 

+   2 

27 

+ 

7 

+   2 

+ 

7 

+   2 

+ 

6 

+   2 

Z 

+ 

44 

+   9 

+ 

44 

+   7 

+ 

44 

+   7 

I  -  III 

+  35 

+  37 

+  37 

de 

.035 

.037 

•  037 

A 

215 

225 

235 

cos 

-  .8192 

-  .7071 

-  .5736 

sin 

-  .5736 

-  .7071 

-  .8192 

ag 

-  .029 

-  .026 

-  .021 

at) 

-  .020 

-  .026 

-   .030 
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Table  XIV     (cont'd) 


I 

III 

I 

III 

I 

III 

6o  - 

70 

2^0-250 

70  - 

80 

250-260 

80  - 

90 

260-270 

38 

+ 

11 

0 

+ 

10 

0 

+ 

10 

+   1 

37 

+ 

13 

0 

+ 

12 

0 

+ 

13 

+   1 

36 

+ 

13 

+  1 

+ 

11 

+  2 

+ 

Ik 

+   3 

35 

+ 

12 

+   5 

+ 

10 

+   5 

+ 

15 

+   6 

L 

+ 

k9 

+   6 

+ 

h3 

+   7 

+ 

52 

+  11 

I  -  III 

+  ^3 

+  36 

+  kl 

d0 

.ok3 

.036 

.(Al 

A 

245 

255 

265 

cos 

-  .1+226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

d| 

-  .018 

-  .009 

-  .00*+ 

dT) 

-  .039 

-  .035 

-  .041 

3^ 

+ 

12 

+   9 

+ 

10 

+   9 

+ 

15 

+   6 

33 

+ 

Ik 

+  10 

+ 

14 

7 

+ 

15 

+   3 

32 

+ 

19 

+   3 

+ 

20 

+   1 

+ 

18 

0 

31 

+ 

20 

0 

+ 

19 

0 

+ 

18 

0 

z 

+ 

65 

+  22 

+ 

63 

+  17 

+ 

66 

+   9 

I  -  III 

+  h3 

+  hi 

+  57 

d0 

.0U3 

.0^7 

.057 

A 

2*4-5 

255 

265 

cos 

-  .4226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

d£ 

-    .018 

-   .012 

-  .005 

dn 

-  .039 

-  .0^5 

-  .057 

30 

+ 

16 

+   1 

+ 

15 

0 

+ 

Ik 

0 

29 

+ 

12 

+   1 

+ 

10 

0 

+ 

10 

0 

28 

+ 

8 

0 

+ 

7 

0 

+ 

7 

0 

27 

+ 

6 

0 

+ 

5 

0 

+ 

5 

0 

E 

+ 

k2 

+   2 

+ 

37 

0 

+ 

36 

0 

I  -  III 

+  to 

+  37 

+  36 

de 

.0*+0 

.037 

.036 

A 

2^5 

255 

265 

cos 

-  .4226 

-  .2588 

-  .0872 

sin 

-  .9063 

-  .9659 

-  .9962 

d| 

-   .017 

-  .010 

-  .003 

dt) 

-   .036 

-  .036 

-  .036 
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Table  XIV     (cont'd) 


II 

IV 

II 

IV 

II 

IV 

90-100 

270-280 

100- 

110 

21 

30-290 

110- 

120 

290-300 

38 

+ 

15 

-      3 

+ 

20 

-      6 

+ 

20 

-       5 

37 

+ 

15 

1 

+ 

20 

-      5 

+ 

20 

-       '+ 

36 

+ 

17 

+       2 

+ 

20 

2 

+ 

20 

0 

35 

+ 

18 

+      6 

+ 

20 

+      5 

+ 

20 

+      5 

L 

+ 

65 

+      k 

+ 

80 

-      8 

+ 

80 

-      1+ 

II  -    IV 

+    61 

+    88 

+     81+ 

d.e 

.06a 

.088 

.081*. 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.1*226 

sin 

- 

.9962 

- 

.9659 

_ 

.9063 

d| 

+ 

.006 

+ 

.023 

+ 

.035 

dT) 

- 

.061 

- 

.085 

- 

.076 

3^ 

+ 

20 

+      5 

+ 

20 

+      5 

+ 

20 

+      5 

33 

+ 

19 

+      3 

+ 

20 

+      1+ 

+ 

20 

+      7 

32 

+ 

18 

+      1 

+ 

18 

+      1 

+ 

18 

+       2 

31 

+ 

16 

0 

+ 

15 

0 

+ 

15 

0 

z 

+ 

73 

+      9 

+ 

73 

+     10 

+ 

73 

+  11*. 

II  -   IV 

+    64 

+  63 

+    59 

d0 

.06k 

.063 

.059 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.  1*226 

sin 

- 

.9962 

- 

.9659 

- 

.9063 

d| 

+ 

.006 

+ 

.016 

+ 

.025 

dT) 

- 

.061*- 

- 

.061 

- 

.051*- 

30 

+ 

13 

0 

+ 

12 

0 

+ 

12 

0 

29 

+ 

10 

0 

+ 

10 

0 

+ 

9 

0 

28 

+ 

7 

0 

+ 

7 

0 

+ 

6 

0 

27 

+ 

5 

0 

+ 

5 

0 

+ 

5 

0 

Z 

+ 

35 

0 

+ 

3* 

0 

+ 

32 

0 

II  -   IV 

+    35 

+    3^ 

+   32 

de 

.035 

.031+ 

.032 

A 

275 

285 

295 

cos 

+ 

.0872 

+ 

.2588 

+ 

.1*226 

sin 

- 

.9962 

- 

.9659 

- 

.9063 

d| 

+ 

.003 

+ 

.009 

+ 

.014 

dt) 

- 

.035 

- 

.033 

- 

.029 
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Table  XIV     (cont'd) 


ii 

IV 

II 

IV 

II 

IV 

120-130 

300-310 

130-140 

310-320 

140-150 

320-330 

38 

+ 

20 

0 

+ 

19 

-   5 

+ 

10 

-   8 

37 

+ 

20 

-   3 

+ 

19 

-   6 

+ 

15 

-   6 

36 

+ 

20 

l 

+ 

19 

2 

+ 

16 

-   3 

35 

+ 

20 

+   2 

+ 

20 

+   1 

+ 

18 

0 

z 

+ 

80 

2 

+ 

77 

-   12 

+ 

59 

-  17 

II  -  IV 

+  82 

+  89 

+  76 

d0 

.082 

.082 

.076 

A 

305 

315 

325 

cos 

+ 

.5736 

+ 

•  7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

d£ 

+ 

.047 

+ 

.063 

+ 

.062 

dT) 

- 

.067 

- 

.063 

- 

.044 

34 

+ 

20 

+   5 

+ 

20 

+   3 

+ 

19 

+   2 

33 

+ 

19 

+   8 

+ 

17 

+   3 

+ 

15 

+   2 

32 

+ 

17 

+   2 

+ 

14 

+   2 

+ 

13 

0 

31 

+ 

Ik 

0 

+ 

12 

0 

+ 

11 

0 

S 

+ 

70 

+  15 

+ 

63 

+   8 

+ 

58 

+   4 

II  -  IV 

+  55 

+  55 

+  54 

&e 

.055 

.055 

.054 

A 

305 

315 

325 

cos 

+ 

.5736 

+ 

.7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

<U 

+ 

.032 

+ 

.039 

+ 

.044 

dT) 

- 

.01+5 

- 

.039 

- 

.031 

30 

+ 

11 

0 

+ 

11 

0 

+ 

10 

0 

29 

+ 

9 

0 

+ 

9 

0 

+ 

8 

0 

28 

+ 

7 

0 

+ 

7 

0 

+ 

6 

0 

27 

+ 

5 

0 

+ 

5 

0 

+ 

5 

0 

Z 

+ 

32 

0 

+ 

32 

0 

+ 

29 

0 

II  -  IV 

+  32 

+  32 

+  29 

d0 

.032 

.032 

.029 

A 

305 

315 

325 

cos 

+ 

.5736 

+ 

.7071 

+ 

.8192 

sin 

- 

.8192 

- 

.7071 

- 

.5736 

d£ 

+ 

.018 

+ 

.023 

+ 

.024 

ai 

- 

.026 

- 

.023 

- 

.017 
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Table  XIV     (cont'fr) 


II 

IV 

II 

IV 

II 

IV 

150- 

l6o 

330-3^0 

l6o- 

170 

3^0-350 

170-180 

350-360 

38 

+ 

k 

2 

+ 

1 

0 

+   5 

+   5 

37 

+ 

7 

-   k 

+ 

3 

0 

+  k 

+   5 

36 

+ 

10 

-   3 

+ 

k 

+   2 

+      1 

+   8 

35 

+ 

11 

1 

+ 

1 

+   2 

0 

+  10 

2 

+ 

32 

-  10 

+ 

9 

+   k 

+  10 

+  28 

II  - 

IV 

+  k2 

+   5 

-  18 

ae 

.OkZ 

.005 

.018 

A 

335 

3^5 

175 

cos 

+ 

.9063 

+ 

.9659 

-  .9962 

sin 

- 

.U226 

- 

.2588 

+  .0872 

ag 

+ 

.038 

+ 

.005 

-  .018 

dl 

- 

.018 

- 

.001 

+  .002 

3^ 

+ 

12 

+   1 

+ 

5 

+   3 

+   1 

+   8 

33 

+ 

15 

+   1 

+ 

10 

+   3 

+   3 

+   6 

32 

+ 

13 

+   1 

+ 

9 

+   2 

+   6 

+   5 

31 

+ 

10 

0 

+ 

8 

+   2 

+   7 

+   5 

£ 

+ 

50 

+   3 

+ 

32 

+  10 

+  17 

+  2k 

II  - 

IV 

+  1^7 

+  22 

-   7 

de 

.0^7 

.022 

.007 

A 

335 

3*+5 

175 

cos 

+ 

.9063 

+ 

.9659 

-  .9962 

sin 

- 

.4226 

- 

.2588 

+  .0872 

d| 

+ 

.0^3 

+ 

.021 

-  .007 

dt] 

- 

.020 

- 

.006 

+   .001 

30 

+ 

9 

0 

+ 

8 

+   2 

+   6 

+      5 

29 

+ 

8 

0 

+ 

7 

+   3 

+   5 

+   5 

28 

+ 

6 

0 

+ 

5 

+   1 

+   5 

+   5 

27 

+ 

5 

0 

+ 

5 

+   1 

+  k 

+   3 

2 

+ 

28 

0 

+ 

25 

+   7 

+    20 

+  18 

II  - 

IV 

+  28 

+  18 

+   2 

de 

.028 

.018 

.002 

A 

335 

3^5 

355 

cos 

+ 

.9063 

+ 

.9659 

+  .9962 

sin 

- 

.4226 

- 

.2588 

-  .0872 

d| 

+ 

.025 

+ 

.017 

+  .002 

dt) 

- 

.012 

- 

.005 

-   .000 
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Table  XV     Summation  of  jtffect*;  o^  Inner  Circle 
and  Zones  27  -  48  for  Laiska 


Zone    Group  Sum  Accumulative  Group  Sum  Accumulative 
"froup     dg"         d£"        &t\"  dT)M 


Outer 

Radius 

(km.) 


Inner 

+ 

0.010 

+  0.010 

-  0.268 

-  0.268 

10.15 

27  -  30 

- 

0.023 

-  0.013 

-  0.383 

-  O.651 

20.09 

31  -  34 

- 

0.062 

-  o\o75 

-  0.601 

-  1.252 

39.67 

35  -  38 

- 

0.072 

-  0.147 

-  0.692 

-  1.944 

77.97 

39  -  42 

+ 

0.161* 

+  0.017 

-  0.667 

-  2.611 

151.9 

43  -  46 

+ 

O.1+68 

+  0.485 

-  0.763 

-  3.374 

291.2 

47  -  48 

+ 

0.060 

+  0.545 

V 

-  0.866 

-  4.240 

399.0 

g 

+ 

0.545 

-  4.240 

100 


kl 

V) 


I 
VD 


v 
V) 


0 
K 

J5 


*  «*J 


S<7/vojj  ?     a-/       fi>oij.r>3~rj3Q 
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7.  V3  Computation  for  Wonderfontein,  Union  of  South  Africa. 

As  mentioned  in  section  7.^1,  both  Bouguer  anomaly  maps  and 

topographic  maps  were  required  for  this  station.  This  situation 

prevailed  from  the  outer  zone,  no.  ^8,  through  zone  no.  31*  The  area 

within  the  circle  described  by  that  zone's  radius,  20.09  km.,  was 

covered  by  a  larger  scale  free -air  anomaly  map  prepared  by  the  writer. 

The  inner  circle  for  this  station  had  a  radius  r  =  5«125  km. 

o 

Figure  l6,  below  is  a  reproduction  of  the  inner  circle.  Using  the 
values  shown,  we  obtains - 

d|"  =  0.0525  (AgN  -  Agg) 

=  0.0525  (99.0  -  88.5)  =  0.0525  (10. 5)  =  +  OV551 

dtl"    =   0.0525    (Agw   -  AgE) 

=  0.0525   (88.2  -  90.2)  =  0.0525   (-2.0)  =-0f.'105 

J  A/ 


99.  o 


-  -S—s^--\  9o. z 


88.S 


Fig.   l6   -  Inner  Circle   -  Wonderfontein. 
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The  effects  of  both  the  inner  circle  and  the  circle-ring  zones 
are  accumulated  in  Table  XVI.  Curves  of  the  effect  on  deflections  of 
increasing  gravity  field  radius  are  shown  in  fig.  17. 

7.^  Comments  on  the  Results  of  Deflection  Computation. 

The  results  shown  in  Tables  XV  and  XVT  and  in  Figures  15  and 
1 7  are  the  computed  effects  of  the  gravity  field  around  each  point 
out  to  a  radius  of  399'0  km.  ^  The  effect  of  the  rest  of  the  world, 
beyond  399' 0  km.  should  not  be  ignored  in  practice.  Time  and  avail- 
ability of  data  do  not  permit  careful  computation  of  this  effect  for 
purposes  of  this  thesis.  However,  an  idea  of  its  magnitude  can  be 
gotten  from  publications  containing  results  of  investigations  into 
this  problem  by  various  authors.  For  instance,  Kaula  (19)>  (20), 
determined  the  RMS  uncertainties  in  the  total  deflection  as  computed 
from  the  Vening  Meinesz  formulas,  first  assuming  perfect  knowledge  of 
gravity  within  given  radii  and  none  beyond,  and  second,  ideal 
distribution  of  given  numbers  of  gravity  stations  out  to  given  radii. 
Taking  Kaula' s  figures  for  the  first  case,  assuming  perfect  knowledge 
within  a  radius  of  200  km.  and  none  beyond,  the  RMS  uncertainty  is 
k"l;     for  500  km.  it  is  3'.'3;  for  1000  km.  it  is  276;  etc.  It  should 
be  noted  that  Kaula' s  computations  were  based  on  information  provided 
by  Hirvonen  (39^- 

Interpolating  for  399  km.  radius,  an  RMS  uncertainty  of  +  3V6 
is  obtained  if  Kaula1 s  information  is  assumed  to  apply.  The  probable 
error  is  then  +  0.67^5  (3V6)  =  +  2?k. 
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Table  XVI  -  Summation  of  Effects  of  Inner  Circle 
and  Zones  23  -  48  for  Wonderfontein. 


Outer 
Zone    Group  Sum  Accumulative  Group  Sum  Accumulative  D  , . 

Group  d£"  d£"  dri"  dT)"  7}^S) 


Inner  +  0.551 

23-26  +  0.515 

27-30  +  0.905 

31-34  +  0.128 

35-38  +  0.175 

39-^2  -  1.884 

43  -  46  -  1.919 

47-48  -  0.75^ 


Z       I    =  -  2.283 


+  0.551 

- 

0.105 

-  0.105 

5.125 

+  1.066 

+ 

0.013 

-  0.092 

10.15 

+  1.971 

+ 

0.212 

+  0.120 

20.09 

+  2.099 

+ 

1.099 

+  1.219 

39.67 

+  2.274 

+ 

i.m 

+  2.667 

77.97 

+  0.390 

+ 

1.139 

+  3.806 

151.9 

-  1.529 

+ 

1.222 

+  5.028 

291.2 

-  2.283 

=  + 

0.493 

+  4.535 

399.0 

*L 

4.535 
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Estimates  of  the  probable  errors  of  the  gravimetric  deflections 
as  computed  from  the  data  given  within  the  field  of  radius  399*0  km. 
must  be  made,  since  perfect  knowledge  of  gravity  within  the  field  can 
not  be  assumed.  The  method  to  be  used  for  this  estimate  is  that  of 
the  U.S.C.  &  G. S.  and  is  described  in  Rice's  paper  on  the  l6  points 
(26).  According  to  this  method,  the  probable  error  in  either  component 
£  or  i),   is  determined  from  the  following  expression: - 

E,2  =  18  Z  (O'.'OOl  E  f 

X>  I* 

where     E    =  Probable  error  in  determining  the  mean  anomaly 

in  any  compartment  in  a  sector  in  mgals. 
E.   =  Probable  error  in  |  or  t;  in  seconds 

In  making  use  of  the  above  equation  for  this  problem,  the 
probable  error  was  estimated  for  each  zone  group  computed.  The 
accumulated  results,  combined  with  the  effect  of  the  inner  circle,  each 
component  of  which  is  assumed  to  be  no  more  certain  than  its  value,  and 
with  the  probable  error  for  the  area  beyond  the  radius  399«  0  km. ,  give 
the  following  (see  tables  XVII  and  XVIII):- 

For  Total  Deflection  9   :- 

Lai  ska       Et  =  ±  2"^2 
Wonderfontein  E.  =  +  2VV7. 

Thus  it  appears  that  for  computations  based  upon  the  gravity 
data  given,  and  within  a  field  of  radius  399*0  km.,  the  probable  errors 
for  the  gravity  field  are  almost  negligible  compared  with  those  re- 
sulting from  the  neglect  of  gravity  beyond  that  radius. 
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Table  XVII  -  Computation  of  Probable  Errors  In  Gravimetric 

Deflections  as  computed  for  the  Gravity- 
Field  within  Radius  399.0  km. 


(1)  (2)  (3) 


E 
Ag 


B 


(4) 


(5) 

E 


(6) 


(7) 


(8) 


0.001  E     (0.001  E   )2  Z(0.001EJ2 
r  r  r 


Zones     °^°s   (mgal)    (n^al)        (cigal) 

|[(3)+00']  0.001(5)        (6): 


for  Group 
(2)    (7) 


LAISKA 
27-30 

31-3^ 
35-38 
39-42 
43-46 
47-48 


4 

4 
4 
4 
4 
2 


1 
1 
1 
5 
5 
5 


0.001 
.001 
.001 
.005 
.005 
.005 


0. 000001 
. 000001 
. 000001 
. 000025 
.000025 
.  000025 


0.000004 
.000004 
. 000004 
. 000100 
. 000100 

. 000050 


W0NDERF0NTEIN 
23-26 


27-30 
31-34 
35-38 
39-42 
43-J46 
47-48 


4 
4 
4 
4 
4 
4 
2 


1 
1 
1 
5 
5 
5 
5 


2 
2 
2 
6 

11 
11 
11 


0.002 

0. 000004 

0. 000016 

.002 

. 000004 

. 000016 

.002 

. 000004 

. 000016 

.006 

.000036 

. 000144 

.011 

. 000121 

. 000484 

.011 

. 000121 

.000484 

.011 

. 000121 

. 000242 

1  100  ft.  @  0.03378  mgal  /  rt. 

2  200  ft.  @  0.03378  mgal  /  ft. 

3  500  ft.  @  0.03378  mgal  /  ft. 


Inner  Circles s- 


Laiska 


Wonderfontein 


i    ■    Kt  -  as 

+  O'.'OIO 

+  0.551 

Et2 

T\       I       E   =  dT) 

-  0.268 

0.000100 

-  0V105 

0.303601 

2 
Et 

0.071824 
0.071924 

0. 011025 

2     "2    2 
6     :   Et  =  dT 

+ 

a.2 

0.314626 
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Table  XVII  (cont'd) 

2  2 

E,   =  18  2(0.001  Er)  ;  E.  =  p.e.  in  either  component 

(|  or  tj)  of  the  deflection 


(9) 

(10) 

(11) 

(12) 

Zones 

18  (8) 

Accumulated 

Acci 

emulated 

(11) 

LAISKA 

* 

27-30 

31-3^ 
35-38 
39-J>-2 
43-46 
47-48 

0.0000Y2 
. 000072 
. 000072 
. 001800 
. 001800 
. 000900 

0.000072 
. 000144 
.000216 
.  002016 
. 003816 
. 004716 

+ 

0V008 
.012 
.015 
.045 
.062 
.068 

Total 
WONDER 

Deflection  6   z 
FONTEIN 

=  2(o.0047l6) 
=  0.009432 
=  0"097 

23-26 
27-30 
31-34 
35-38 
39-^2 
43-J46 
47-48 

0.000288 
.  000288 
. 000288 
.002592 
. 008712 
.008712 
. 004356 

0.000288 
. 000576 
. 000864 
. 003456 
. 012168 
. 020880 
. 025236 

+ 

o'ioi7 

.024 

.029 

.059 
.110 
.145 
.159 

Total  : 

Deflection  Q   1 

=  2(0.025236) 
=  0.050472 
=  0V225 
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Ta,ble  XVIII 


Probable  Errors  in  Gravimetric  Deflections  9 


for  Combined  Zones 


Zones 


E. 


Accumulated 
o 

E  * 


Accumulated 

E^ 


LAISKA 


Inner  Circle    (r     =  10.25  km. ) 

0.0719 

0. 0719 

j'.'268 

27 -U8  (10.25  -  399-0  km.) 

0.009^ 

0.0813 

0':285 

Beyond  399-0  km.    !    (2'.'J+)2 

5.76 

5.8IH3 

i*.'lf2 

W0NDERF0NTEIN 


Inner  Circle  (r  =  5«125  km.) 
23-26  (5.125  -  399-0  km.) 
Eeyond  399-0  km.  :  (2*.*U)2 


0.31^ 

0.31^ 

07561 

0.0505 

0.3651 

o':6ok 

5.76 

6.1251 

27  Vf 
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T«5  Approximate  Correction  to  Finland  -  South  Africa  Arc  from 

Gravimetric  Deflections 

Making  use  of  the  values  ©btained  for  the  gravimetric  deflect- 
ions of  Laiska  and  Wonderfontein  in  section  7 ok,   disregarding  the  un- 
certainty arising  from  the  fact  that  computation  was  not  carried  out 
"beyond  399  km0 ,  a  correction  to  the  length  of  the  meridian  arc  between 
the  points  can  be  made.  From  this,  a  rough  idea  of  the  corresponding 
correction  to  the  equatorial  radius  can  be  learned. 

Geodetic  and  astronomic  data  for  the  two  points  in  question 
was  obtained  from  the  Army  Map  Service,  and  was  part  of  the  data  used 
in  that  agency1 s  computation  of  "A  New  Determination  of  the  Figure  of 
the  Earth  from  Arcs"  (See  Chovitz  and  Fischer  (k)).     Unfortunately  the 
data  carries  a  security  classification  for  values  having  accuracy 
better  than  one  minute  of  arc.  Because  of  this,  the  actual  computation 
figures  cannot  be  published  here.  However,  the  results  given  reflect 
the  data  furnished  to  the  writer  by  AMS  without  violating  security 
regulations. 

Using  latitude  functions  tables  (either  (3*0  or  (36))  for  the 
International  Ellipsoid,  upon  which  the  given  data  was  based,  the 
length  of  the  meridian  arc  between  Laiska  and  Wonderfontein  for  the 
geodetic  latitudes  and  for  the  latitudes  corrected  by  the  free -air 
gravimetric  and  astro-geodetic  deflections,  were  computed  (See  figure 
18).   In  effect,  the  two  figures  thus  obtained  were  (l)  the  measured 
length  of  the  arc,  based  on  the  International  Ellipsoid,  and  (2)  the 
corrected  length  of  the  arc  based  on  astro-gravimetric  data.  The 
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following  figures  give  pertinent  details :- 

Station  Latitude         Meridian  Arc 

Laiska  cp  =  6k°   03'  (  )"  N       M 

Wonderfontein  cp  *  25°  48*  (  )"  S       M 


Acp  =  89°  51'  (  )" 


(Mj.  +  Mw)  =  9,961, 


kkQ 


m. 


Laiska  cp  =  6k°   03*  (  )"   N      (M_) 

Wonderfontein     cp  =  25°  48'  (  )"   S      (M  ) 

g         v  '  x  w  * 


Acpff  =  89°  51«  (  )M 

(^  +  Mw)g  =  9,961,534  m. 

where  ]VL  and  M  refer  to  the  meridian  distance  from  the  equator  to  the 

point  in  question,  cp  is  the  geodetic  latitude,  and  cp  the  gravimetric - 

S 

ally  corrected  latitude.  Thus  cp  =  cp  +  (cp'  -  cp)  -  £  . 

o  o 

From  the  above  figures,  we  see  that  the  total  correction  for 
the  arc  is  +  86  meters.  Assuming  the  arc  to  be  a  meridian  arc  of  a 
sphere  the  size  of  the  earth,  we  can  thus  obtain  a  rough  correction  to 
the  sphere ' s  radius : - 

6R  =  £-  x  86 

&p 

57°o 
=  rf-fy-  x  86  =  +  55  meters  ^  8a 

89^85 
Admittedly,  this  is  only  a  rough  approximation,  but  is  close 
enough  to  the  exact  value  for  the  data  used,  for  illustrative  purposes. 
Besides  this,  the  accuracy  of  the  gravimetric  deflections  is  of  doubt- 
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ful  value,  because  of  the  limited  gravity  field  used. 

Any  person  who  is  lucky  enough  to  he  in  possession  of  data 
such  as  is  used  here,  hut  for  many  points  along  an  arc  of  this  length 
would  proceed  along  the  lines  of  section  7.2,  and  would  doubtless 
obtain  results  of  some  significance. 


LAI  SKA 


MER.IOIAN 


WONOERFONTZlN 


Fig.   18 
The  Finland  -  South  Africa  Arc 
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Chapter  8  -  Closing  Remarks 

It  would  be  presumptuous  on  the  part  of  the  writer  to  title 
this  chapter  "Conclusions",  because  there  can  hardly  be  any  conclusion 
to  something  that  has  not,  in  a  manner  of  speaking,  even  begun.  The 
lack  of  readily  available  data  for  pursuing  a  project  of  such 
magnitude  as  determining  the  best  dimensions  for  the  earth  ellipsoid, 
prevents  a  significant  contribution  to  the  science,  even  if  the  ability 
to  put  such  data  to  proper  use  were  present. 

What  conclusions  can  be  made,  must  be  based  on  what  was  al- 
ready known  beforehand,  and  indeed  could  not  be  made  without  such 
knowledge,  because  of  the  nature  and  amount  of  the  data  used.  The 
first  two  problems  solved  in  Chapter  7  are  what  are  referred  to  here. 
The  use  of  only  astro-geodetic  deflections  of  the  vertical  is  admitted- 
ly an  incomplete  beginning  to  a  solution,  so  could  hardly  be  expected 
to  yield  as  good  results  as  with  the  use  of  the  gravimetric  deflections. 
It  is  unfortunate  that  gravimetric  deflections  were  not  available  for 
all  the  points  used  for  the  arc  of  the  98th  meridian.  A  direct  com- 
parison of  methods  could  then  have  been  made.  It  would  have  been 
almost  as  meaningless  to  have  computed  the  size  of  the  earth  from  only 
the  astro-geodetic  deflections  of  the  l6  points  in  the  South  Central 
U.S.,  because  of  the  small  area  covered.  As  it  was,  the  computation 
by  the  astro-gravimetric  method  in  the  case  of  the  16  points,  probably 
owed  much  of  its  "success"  to  the  fact  that  the  data  given,  as 
prepared  by  Mr.  Rice,  was  carefully  computed  by  him  and  fairly 
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accurate. 

It  is  hoped  that  the  material  presented  in  this  thesis  might 
provide,  by  its  systematic  presentation  of  the  procedure  to  be 
followed  in  problems  of  this  nature,  the  means  whereby  any  interested 
reader  might  pursue  the  problem  further,  especially  when  in  the  future 
more  gravimetric  deflections  of  the  vertical  can  be  expected  to  be 
available.  The  advent  of  the  high  speed  electronic  computer  in 
geodetic  use  will  hasten  that  day.  The  manual  method  of  computing 
gravimetric  deflections  is  so  laborious  that  one  person,  in  the  work 
of  preparation  of  a  thesis,  could  not  be  expected  to  produce  the 
required  computations  in  adequate  numbers  and  for  those  stations  in 
locations  favorable  to  solution  of  the  problem  of  determining  the  size 
of  the  earth  ellipsoid.   It  is  believed  that  the  example  deflection 
computations  included  in  this  thesis  demonstrate  this  clearly. 

It  is  regretted  that  more  than  mere  mention  of  two  important 
methods  used  for  solution  of  the  problem  at  hand  could  not  have  been 
made  in  this  thesis.  One  is  the  determination  of  the  figure  of  the 
earth  from  geoid  heights,  and  the  other  is  the  determination  of  the 
flattening  from  observations  of  the  artificial  earth  satellites. 
Time,  more  than  space,  prevented  their  inclusion,  especially  since  the 
writer  tends  to  be  somewhat  verbose  in  his  explanations. 
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